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This publication describes the mathematical and service subpro- 
grams contained in the library supplied with Basic FORTRAN Iv and 
FORTRAN Iv compilers. Detailed information on each mathematical 
routine is provided: its algorithm, speed, accuracy, range, and 
error conditions. Information on the use of library subprograms 
in either a FORTRAN or assembler language program is also provided. 


~ Preface 


The purpose of this publication is to describe the 
mathematical and service subprograms in the FORTRAN’ 


1v library supplied with Basic FORTRAN Iv (os) and 

FORTRAN Iv (os), As an aid to the programmer in his 

use of this publication, the contents of each chapter is 

described briefly below: 

1. “Introduction” describes the three types of subpro- 
grams in the FoRTRAN library (FORTLIB)and de- 
fines their use in either a FORTRAN or an assembler 
language program. 

2. “Mathematical Subprograms” describes the subpro- 
grams which perform computations frequently 
necded by the programmer. A mathematical sub- 
program is invoked explicitly whenever one of its 
entry names appears ‘in a- source statement, or im- 


plicitly through use of certain aotation in the source 


statement. 

3. “Service Subprograms” contains infesuaeon about 
those subprograms which perform. utility functions 
or test machine indicators. 

4, “Algorithms” contains information about the method 
used in the library to compute a mathematical func- 
tion and describes the effect of an argument error 
upon the accuracy of the answer returned. 


5. “Performance Statisticn” gives agmMacy ¢ end tiraing 


statistics for the expiicn calied ana dictaatical sub: 
programs. 


SECOND EDITION (September 1972) 


6. Theappendixes provide a list of diagnostic messages, 
a list of module names, a sample storage printout, 
storage estimates, and information for the assembler 
language programmer. 

It is assumed that the reader is familiar with one of 
the following publications: 
IBM System/360 F ORTRAN IV Language, Form 

C28-6515 . 
IBM System/360 Basic FORT RAN IV Language, 

Form C-28-6629 
IBM System/360 Operating System: Assembler Lan- 

guage, Form C28-6514 
In addition, references are made within this publica- 

tion to information contained in ne following publica- 
tions: 
IBM System/360 Principles of Operation, Form 

A22-6821 
IBM/360 Operating System: Supervisor and Data 

Management Macro-Instructions, Form C28-6647 

IBM System/360 Operating System: Basic FORTRAN 

IV (E) Programmer’s Guide, Form C28-6603 
IBM System/360 Operating System: FORTRAN IV 

(Gand H) Programmer’s Guide, Form C28-6817 
Standard mathematical notation is used in this 

publication. The reader is expected to be familiar 

w:th_ this noiation and y es commen mathematical 


OT inology. 


This is a reprint of GC28-6818-0 incorporating changes released in Technical Newsletter GN28-0589, 


dated Jnne 1, 1970. 


This edition corresponds to Release 19 of the IBM System/360 Operating System. 
Changes are periodically made to the specifications herein; any such changes will be reported in 


subsequent revisions or Technical Newsletters. 


Requests for copies of IBM publications should be made to your IBM representative or to the IBM 


Branch Office serving your locality. 


A form for readers’ comments is provided at the back of this publication. If the form has been 
removed, comments may de adére- +d to IBM Corporation, Programming Publications, 


1271 Avenue of the Americas, New York, New York 10020. 
© Copyright International Business Machines Corporation 1969 
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The rortaax a library contains three types of sub- 
programs: (1) mathematical functions such as sin and 
sort, (2) service subprograms such as pump and EXIT, 
and (3) interface routines such as 1con and BERR. 

The library subprograms may be used in either a 
FORTRAN or an assembler language program. (Appen- 
dix A contains calling information for the assembler 
language programmer.) In Fortran. calls to the library 
subprograms are cither at the programmer's request 
or in response to the program requirements. Subpro- 
grams required by the program being compiled are 
provided by a linkage editor or loader. which takes 
the subprograms from the library. 


Mathematical Function Subprograms 


Many commonly used mathematical functions or 
calculations are provided by the ForTRAN Iv language. 
Programmer requests for these calculations are satisfied 
by a compiler in one of two ways: by inserting machine- 
language code directly within the object module; or 
“by inserting an extemal reference to a library subpro- 
gram that contains the necessary code. and having the 
linkage editor include that subprogram in the load 
module. The first method is called in-line, the second 
out-of-line. 

This publication discusses only out-of-line routines— 
these mathematical functions contained in the library. 
The American National Standards Institute (anst) de- 
fines several arithmetic functions—such as absolute 
value (ans), positive difference (pive). and transfer of 
sign (siGN )—as intrinsic functions. For the most pari. 


Introduction 


code for these functions is inserted in-line by the 
FORTRAN compiler at the point in the source module 
where the function's symbolic name is used. Conse- 
quently, they are not discussed in this book. However, 
the following ansi-defined intrinsic functions have been 
implemented as a part of the FoRTRAN IV library and 
are provided out-of-line for all systems: .fax0/ INO, 
ANEAXO/ AMINO, MANI/MINI, AMANI/AMINI, and DMAXx1/ 
parini. These are, therefore, documented in this publi- 
cation. In Basic FORTRAN Iv only, four other intrinsic 
functions (\fop, AMOD, AINT. and IFIX/INT/IDINT ) are 
provided out-of-line as a part of the library and are also 
covered here, 


Service Subroutine Subprograms 

Each of the service subprograms corresponds to a sub- 
routine form as defined by a sUBROUTINE statement in 
FORTRAN source module. These subprograms perform 
machine indicator tests and, utility functions and may 
or may not return a value to the calling module. These 
subroutines are discussed in this publication. 


Interface Routines 


The library contains certain input/output and error 
processing routines that act as interfaces with the 
compiled program and the operating system. Fre- 
quently. the mathematical and service subprograms 
require assistance from these routines for input/output, 
interruption. and error processing. Storage estimates 
for these routines are included in Appendix B. 
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Mathematical Subprograms 


The mathematical subprograms supplied in the = power) appears within a source language statement 
FORTRAN library perform computations frequently (see Table 6). 
needed by the programmer. The mathematical subpro- The following text describes the individual mathe- 
grams are called in two ways: explicitly, when the — matical subprograms and explains their use in a 
programmer includes the appropriate entry name ina = FORTRAN: program. Detailed information about the 
source language statement (see Table 1); andimplicitly, | actual method of computation used in each. subpro- 
when certain notation (e.g.. raising a number to a — gram, the’ performance of the subprogram. interrup- 
tion and error procedures. and storage estimates can 
be found elsewhere in this publication. 


Table 1. Explicitly Called Mathematical Subprograms 


General Function Specific Function Entry Name(s) 


Logarithnic and exponential subprograms Exponential 
(described in Table 2) 

























ALOG, ALOGIO 
DLOG, DLOG10 

CLOG* 
CDLOG* 


SORT 


Logarithmic. common and natural 
DSQRT 
CSQRT® 


Square root 
CDSQRT* 


Aresine and arccosine ARSIN*, ARCOS* 
DARSIN*, DARCOS* 

Arctangent _ | ATAN, ATANQ® 7 
DATAN, DATAN2® 


Sine and cosine SIN, COS 
DSIN, DCOS 
CSIN*, CCOS® 
CDSIN*,. CDCOS* 


Tangent and cotangent TAN*, COTAN*® 
DTAN*, DCOTAN* 
















Trigonometric subprograms 
(described in Table 3) 























































Hyperbolic function subprograms Hyperbolic sine and cosine SINH*, COSH* 
(described in Table 4) DSINH*, DCOSH* 








Hyperbolic tangent TANH 


DTANH 
Absolute value CABS* 
CDABS* 
ERF*, ERFC* 
DERF*, DERFC* 
Gamma and log-gamma GAMMA*, ALGAMA*® 
DGAMMA‘*, DLGAMA* 


Maximum and minimum value AMAXO, AMINO, MAXO, MINO 
AMAXI1, AMINI, MAX], MINI] 
Modular arithmetic MOD 
AMOD, DMOD 


DMAX1, DMIN1 
Truncation AINT 
; INT, IDINT 


*Not available in Basic FORTRAN IV (OS) 






Miscellaneous subprograms 
(described in Table 5) 









Explicitly Called Subprograms 

Each explicitly called subprogram performs one or 
more mathematical functions. Each mathematical fune- 
tion is identified by a unique entry name. 

A subprogram is called whenever the appropriate 
entry name is included in a FORTRAN arithmetic expres- 
sion. The programmer must also supply one or more 
arguments. These arguments follow the entry name 
and are separated by commas; the list of arguments is 
enclosed in parentheses. 

For example, the source statement: 

RESULT = SIN (RADIAN) 
causes the sine and cosine subprogram to be called. 
The sine of the value in RADIAN is computed and the 
function value is stored in RESULT. 

In the following example, the square root subpro- 


Table 2. Logarithmic and Exponential Subprograms (Part I of 2) 


Entry Sub- 
Name set 


General 


Function Definition 


Common 
and natural 
logarithm 


y=log.x or 
y=Inx 


ne 


ne 
% 


me a) 
Sho ‘ 


y=PV log. (7) 
See Note 2 


y=PV log. (z) 
bee Note 2 


es | y=e 





B 
“< 
a 
% 
“ 


“< 
oi 
2% 
So 
® 
rx) 


< 
% 


_ = 


Notes: (See end of table. ) 


cela s) 


SOMPLEX *§8 


MPLEX °16 [xi SS 174.673 
eee < (250° #) 


gram is called to compute the square root of the value 
in aMINT. The function value is then added to the value 
in sTOCK and the result is stored in ANs. 
ANS = STOCK + SQRT (AMNT) 
The explicity called subprograms are described in 
Tables 2 through 5. The following information is pro- 
vided: 


zeneral Function: This column states the nature of 
the computation performed by the subprogram. 

Entry Name: This column gives the entry name that 
the programmer must use to call the subprogram. A 
subprogram may have more than one entry name, the 
particular entry name used depends upon the compu- 
tation to be performed. For example, the sine and 
cosine subprogram has two entry names: SIN and Cos. 
If the sine is to be computed, entry name sin is used; 
if the cosine is to be computed, entry name cos is used. 


Function Value 
Type’ and 
Range‘ 


REAL *4 
y 2 -180.218 
y S 174.673 


REAL *8 
y = -—180.218 
y S 174.673 


@ 


iv 
Ge 


COMPLEX *8 

y: = —180.218 
y: & 175.021 
-sSy:Sr 


COMPLEX *16 
vy: 2 — 180.218 
¥: & 175.021 
—wFT s yz = T 


REAL *4 
OSySy7 


REAL *§ 
OSysv7 


COMPLEX *8 
—-YSy.y:57 


COMPLEX *16 
-7s yi, Ys 7 


x & 174.673 
x S 174.673 


x: S 174.673 
att < (28° w) 
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Table 2, Logarithmic and Exponential Subprograms (Part 2 of 2) 


Function Value 
General Eat Type! and Error 
Function Name Definition : Ranget Code 


Square root 


y= Vxor (2 REAL *8 
y=x'? . Osysy""" 


y= Vzor any COMPLEX iia *8 

y=z'? argument 0s y: S 1.0987 (4'*) 

See Note 3 tal = 1.0987 (7'7) 

y=Vzor 1 COMPLEX *16 | any COMPLEX CONPLES *16 

y=z'? argument = = vy, S 1.0987 (y'?) 

See Note 3 vz) S 1.0987 (7'*) 
Notes: ~ 


* REAL*4 and REAL*8 arguments correspond to REAL and DOUBLE PRECISION arguments, respectively, in Basic FORTRAN IV. Com- 
plex arguments cannot be used in a Basic FORTRAN IV program. 


* py = principal value. The answer given (y, + y:') is that one whose imaginary part (y:) lies between — x and +. More 
specifically: -— x < y: S 7m, unless x, << O and x: = —0, in which case, y: = — 1. 


= is a complex number of the form x, + x: i. 
y = 16% (1 — 16°) for regular precision routines, and 16% (1 — 16°''). for double-precision routines. 
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Subset: This column indicates those subprograms 
that belong to the Basic rorTnan tv library. Unless 
otherwise indicated, all such subprograms also belong 
to the FORTRAN IV libraries. 


Definition: This column gives a mathematical equa- 
tion that represents the computation. An alternate 
equation is given in those cases where there is another 
way of representing the computation in mathematical 
notation. (For example, the square root can be repre- 
sented either as y = Vx or y = x'/2,) 


Argument Number: This column gives the number 
of arguments that the programmer must supply. 


Argument Type: This column describes the type and 
length of the argument(s). INTEGER, REAL, and COMPLEX 
represent the type; the notations *4, *8, and *16 repre- 
sent the size of the argument in number of storage 
locations. 

Argument Range: This column gives the valid range 
for arguments. If an argument is not within this range, 
an error message is issued and execution of the load 


Table 3. Trigonometric Subprograms. (Part I of 2) 


ARSIN 


aia a oa 






Arcsine and 
arccosine 








Arctangent 


cosine 


Notes: (See end of table. ) 





> 


REAL °4 jx] < (2!8 ¢ x) 
(in radians ) 


(in radians ) 
*8 
(in radians ) 


G 


module is terminated unless the extended error han- 
dling facility was specified at system generation (see 
FORTRAN IV (G and H) Programmer's Guide listed 
in the Preface for a full description of this facility). 
Appendix C of this publication contains a description 
of the error messages. 


' Function Value Type and Range: This column de- 
scribes the type and range of the function value re- 
turned by the subprogram. Type notation used is the 
same as that used for the argument type. Range symbol 
y = 16*3(1—16-*) for regular precision routines, and 
16%3(1—16-'4) for double-precision routines. 

Error Code: This column gives the number of the 
message issued when an error occurs. Appendix C 
contains a description of the error messages. 

Throughout this publication, the following approxi- 
mate values are represented by 2?° ¢ a and 25° « 7: 

218 « gr = 82354966406249996D + 06 

250 « a = .35371188737802239D + 16 

Detailed information for the assembler language 
programmer is given in Appendix A. 


Function Value 
t 
Argument(s) Sepp eere anon Type’ and Error 
Range* Code 


REAL °4 (in radians) | 257 


REAL *4 (in radians) 
OSysr 


REAL *8 (in radians) | 267 

nee yy ae 
go 

REAL °8 (in radians) 

OSysr 

REAL °4 (in radians) 


~~<ys> 


2 


REAL *4 (in radians) | 
-*3r<ysr 


REAL °8 (in radians) 


e 


E 


any REAL argument 


any REAL arguments 
(except 0, 0) 


any REAL argument 


any REAL arguments 
(except 0, 0) 







REAL *8 (in radians) 









ix] < (25° x) 
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7 ATeimenia) Function Value 
General Entry Sub- Definition eee ee Type’ and Error 
Function Name set ‘| No. ‘Type! Range‘ Code 
Sine and CSIN y=sin (z) COMPLEX *8 { [x] < (278 © w) COMPLEX *8 274, 
cosine See Note 2 (in radians) | |x:| S 174.673 - -ySy,yS%7 275 
(continued ) 


an — ees mee, ed 


CCOSs 


y=cos (z) COMPLEX *8 
See Note 2 (in radians ) 


Xr 


Ixs| < (28 w) COMPLEX *8 
i 
Xs 


<= 174.673 -ySyny2357 





ote CORRE ee me 





< (25% wr) . COMPLEX *16 
S 174.673 —ySyny547 


y=sin (z) COMPLEX * 16} |x: 
See Note 2 (in radians ) | |x: 


erry cow reee eee we, 


CDSIN 


ww 
~] 


























CDCOS y=cos (z) COMPLEX *16]|x:] < (25° x) COMPLEX *16 284, 
See Note 2 (in radians) }|x2| S 174.673 —-ySy,yS7 285 
Tangent TAN y =tan (x) 4 REAL *4 
and (in radians ) | See Note 3 
cotangent 
COTAN y=cotan (x) 1 |reactd — [fx] < (2! x) REAL *4 258, 
(in radians) {See Note 3 -ySyS7 259 
DTAN y =tan (x) 1 REAL *8 Ix| < (28° « x) REAL *§ 268, 
(in radians) |See Note 3 —~ySyfyr 269 
DCOTAN |No_ | y=cotan (x) 1 REAL *8 [x] < (25% w) REAL °8 268, 
(in radians) |See Note 3 —“Y¥SyYSy 289 


rns a aw: 


NOrEs: 


1 REAL *4 and REAL °8 arguments correspond tn REAL and DOUBLE PRECISION arguments, respectively, in Basic FORTRAN 1v. Com- 
plex arguments cannot be used in a Basic FORTRAN Iv program. 
* z is a complex number of the form x: + xsi. 


* The argument for the cotangent functions may not approach a multiple of +; the argument for the tangent functions may not 
approach an odd multiple of «/2. 


“y = 16% (1 — 16%) for regular precision routines, and 16° (1 -- 16*) for double-precision routines. . ace 


Peete te TOS 








Table 4, Hyperbolic Function Subprograms 










Ee 














Function Value 


General Entry revert ) Type’ and 
Function Name. : . Range" 


Hyperbolic 
sine and 







Norss: 
: REAL 4 and REAL *8 arguments correspond to REAL and DOUBLE PRECISION arguments, respectively, in Basic FORTRAN Iv. Com- 
plex arguments cannot be used in a Basic FORTRAN Iv program. 


* + = 16" (1 — 16°) for regular precision routines, and 16% (1 — 16°") for double-precision routines, 
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Table 5. Miscellaneous Mathematical Subprograms (Part I of 2) 


General | Entry Sub. 


Function Value 
Type’ and 
Function Name set Definition No. Range* 


Absolute COMPLEX *8 ] any COMPLEX 


value argument 
COMPLEX *16 



























See Note 2 


y=|z{= (x? + x7)" any COMPLEX 
argument 


See Note 2 






Error 
function 







any REAL 
argument 















any REAL 
argument 





Gamma and 
log-gamma 





x > 2-752 and 
x < 57.5744 
















x > O and 
x < 4.2913 « 1078 


REAL *4 


y=log. F (x) or 
? —0.12149 =ySy7 





yolos. f u™! e* du 




















x > 27°52 and 
x < 57.5744 





DGAMMA No 
o 








DLGAMA N 








y=log. (x) or x > Oand 
eo 


x < 4.2913 © 1078 


REAL *8 
—0.12149 =y=y7 










y=bog. f u™te“ du 






Maximum 
and 

minimum 
values 


y=max (xr,... 


















y=min (x;,... 
y=max (x1,... 
y=min (xi,... 


y=max({x,... 


Notes: (See end of table. ) . 
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Table 5. Misce:slaneous Mathematical Subprograms (Part 2 of 2) 


ion Val 
Géneral Entry Sub- Argument(s) Too Error 
Function Name set Definition | No. | Typet | Range | ~ Range* Code 
MINI Yes 


y=min (x1,...,Xn) 22 | reat *4 any REAL INTEGER *4 
arguments 
AMAX1 =max (m,... 


2 
2 














Maximum 
and 
minimum 
values 
(continued ) 





nd 
& 


REAL *4 





any REAL 















IV 


i 


y=min (x:,..., Xa) *4 





nd 
SF 


Yes ly=max (x,... 


any REAL 
ar 








y 
guments 


REAL *8 any REAL 
arguments 
X2 0 
See Note 5 

REAL *4 X2 0 
See Note 5 

REAL *8 x2 % 0 REAL *8 
See Note 5 

REAL *4 ‘| any REAL 

argument 


Yes 


” 
Im 





Modular 
arithmetic 


y =x: (modulo xe) 
See Note 4 


y=x: (modulo x:) 
Note |SeeNote4 | 
3 


Note 


7) 
4 


> “Ss > 
3 |: 
S iS = 





y=x; (modulo xz) 
ee Note 4 


Zw 
os 8 





Truncation |AINT y = (signx) en 
Note |where n is the largest 


integer S |x| 


n 
$ 












o 











INT See ly = (signx) ¢n REAL *4 
Note jwhere n is the largest 
3 integer S |x| 
IDINT See ly = (signx)¢n REAL *8 





where n is the largest 


Note 
3 finteger S |x] 





Norges: 


" REAL®4 and REAL®8 arguments correspond to REAL and DOUBLE PRECISION arguments, respectively, in Basic FORTRAN Iv. Com- 
plex arguments cannot be used in a Basic FORTRAN IV program. 


* Floating-point overflow can occur. 


* The coding that performs this function is out-of-line in Basic FORTRAN Iv (os) and in-line in FORTRAN Iv. Out-of-line coding 
is taken from the ForTnan library by the linkage editor or loader and processed with the calling module. In-line coding is 
inserted by the rorTRAN compiler at the point in the source module. where the function is referenced. This means that the 
in-line functions are invoked in ForTRAN Iv by using the appropriate entry name, but that they are not part of the library. 
In-line functions are described in the ForTRAN 1v language publications listed in the Preface. 





* The expression x: (modulo xz) is defined as x1 — [ x | ° xz, where the brackets indicate that an integer is used. The largest 


Xz 
integer whose magnitude does not exceed the magnitude of = is used. The sign of the integer is the same as the sign of = : 


* if Xs = 0, then the modulus function is mathematically undefined. In addition, a divide exception is recognized and an inter- 
ruption occurs..(A detailed description of the interruption procedure is given in Appendix C.) 


* y = 16" (1 — 16*) for regular precision routines, and 16° (1 — 16") for double precision routines. 
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implicitly Called Subprograms 


The implicitly called subprograms are executed as a 
result of certain notation appearing in a FORTRAN 
source statement. When a number is to be raised to a 
power or when multiplication or division of complex 
numbers is to be performed, the rorrRAN compiler gen- 
erates the instructions necessary to call the appropriate 
subprogram. For example, if the following source 
statement appears in a source module, 


ANS = BASE**EXPON 


where BASE and EXPON are REAL*4 variables, the 
FORTRAN compiler generates a reference to FRXPR#, the 


Table 6. Implicitly Called Mathematical Subprograms 





















Raise an integer 
to an integer 
power 















Raise a real FRXPI# 
number to an 


integer power 











vo 


FDXPI# 








Raise a real FRXPR# 
number to a 


real power 





FDXPD# 












Raise a complex 
number to an 
integer power 


FCDXI# 





FCXPI# 


NotEs: 

1. This name must be used in an assembler lan 
must be included. 

2. 

3. 


arguments cannot be used in a Basic FORTRAN IV program. 
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licit Function 
General Entry' tapes Argument(s) Value ee 
Function Name Reference’ Type® € 
wey and y = 21° 2s COMPLEX °16 
divide complex 
numbers ee eee eee : 
y = 21/2: COMYILEX *16 COMPLEX *16 


, 
EES EE 
s 
s 
ry 
eee 
Cet + eee cate amet [ 


guage program to call the subprogram; the character # is a part of the name and 


This is only a representation of a FoRTRAN statement; it is not the only way the subprogram may be called. 
REAL °4 and REAL °8 arguments correspond to REAL and DOUBLE PRECISION arguments, respectively, in Basic ronTRAN Iv. Complex 


entry name for a subprogram that raises a real number 
to a real power. . 

The implicitly called subprograms in the FORTRAN 
library are described in Table 6. The column headed 
“Implicit Function Reference” gives a representation 
of a source statement that might appear in a FORTRAN 
source module and cause the subprogram to be called. 
The rest of the column headings in Table 6 have the 
same meaning as those used with the explicitly called 
subprograms. Algorithms for implicitly called subpro- 
grams are given in the chapter “Algorithms.” Additional 
information for assembler language programmers is 
given in Appendix A. 































COMPLEX *8 


COMPLEX *8 


i = INTEGER °4 
j 7 INTEGER *4 






a = REAL °4 
j = INTEGER *4 


a = REAL *§ 





Z = COMPLEX *16 COMPLEX °16 
j = INTEGER °4 
Z = COMPLEX *8 COMPLEX *8 


j = INTEGER *4 


For subprograms that involve exponentiation, the 
action taken within a°subprogram depends upon the 
types of the base and exponent used. Tables 7 through 
10 show the result of an exponentiation performed with 
the different combinations and values of base and ex- 
ponent. In these tables, I and J are integers; A and B 
are real numbers; C is a complex number. 


Table 7. Exponentiation with Integer Base and Exponent 
















Exponent (J) 


I 
need EC 
I>1  |Compute the | Function Function 
function valug value = 1 value = 0 
Compute the | Function Function 
function valuel value = I value = ] 
Function Error message | Error message 
value = 0 241 241 
number, function 
value = 1. 
I< —1 [Compute the | Function Function 
function value] value = 1 value = 0 
















Compute the | Function If J is an odd 
function value] value = 1  —_‘} number, function 
value = —1. 


If J is an even 












Exponent (J) 


[aso [eo] i<0 
A>0O _ |{Compute the Function Compute the 
. function value | value = 1 function value 
A=0O {Function Error message | Error message 
value = 0 _ 242 or 242 or 
243 243 
A<0 [Compute the Function Compute the 
function value | value = 1 function value 





Foy! 
AA 
© 

fo 







Tahle 9. Exponentiation with Real Base and Exponent 


Exponent (B) 


pee Bo [ Beo 


A>0O _ |Compute the Function | Compute the 
function value | value = i function value 


Function Error message | Error message 
value = 0 244 or . 2A4or 


245 245 


Error message | Function Error message 
253 or value = 1 253 or 
263 - 263 


Table 10. Exponentiation with Complex Base and Integer 
Exponent 


C=P+Qi 
> 0 and 


Or 
wy 
VV 2 
oe 
9 


oo!|o 
Qu 


and 


Ovlovwsloe 
VIilIAVIayVv 
Qopjfeosc 
Bi 8 
[ow 


Orv 

oo 
> 
3 
o¥) 


wlowvlor 
5 
tov) 


_@) 

NAILVATA Ul 
SlLOoOOolToeooloo 
J & 
=) =} 

Qs a, 


nd 








Exponent (J) 





Compute the | Function Compute the 
function value} value=1 + 0; | function value 
Compute the | Function Compute the 
function value | value=1 + 0i | function value 
Compute the | Function Compute the 
function value | value= 1 + 0 | function value 


Compute the | Function Compute the 


Ex r<0 


function value | value=1 + Oj | function value 


Function 

value= 1 + 0i | function value 

Function Compute the 
function value | value=1 + 0 [function value 









Compute the 
function value 


Compute the 





Compute the {Function Compute the 
function value | value=1 + 0i [function value 


Compute the | Function Compute the 





Mathematica] Subprograms 


function value | value=1 + Oi [function value 
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Service Subprograms 


The service subprograms supplied in the FORTRAN 
library are divided into two groups: one group tests 
machine indicators and the other group performs 
utility functions. Service subprograms are called by 
using the appropriate entry name in a FORTRAN lan- 
guage CALL statement. 


Machine Indicator Test Subprograms 


The machine indicator subprograms test the status of 
pseudo indicators and may return a value to the calling 
program. When the indicator is zero, it is off; when the 
indicator is other than zero, it is on. In the following 
descriptions of the subprograms, i represents an integer 
expression and j represents an integer variable. 


Pseudo Sense Light Subprogram 
Entry Names: SLITE/SLITET 
This subprogram is used to alter, test, and/or record 
the status of pseudo sense lights. Either of two entry 
names (SLITE Or SLITET ) is used to call the subprogram. 
The particular entry name used in the CALL statement 
depends upon the operation to be performed. 

If the four sense lights are to be turned off or one 
sense light is to be turned on, entry name sLITE is used. 
The source language statement is: 


CALL SLITE (i) 


where i has a value of 0, 1, 2, 3, or 4. 

If the value of i is 0, the four sense lights are turned 
off; if the value of i is 1, 2, 3, or 4, the corresponding 
sense light is turned on. If the value of i is not 0, 1, 2, 
3, or 4, error message 216 is issued and execution of 
this module or phase is terminated. (This error han- 
dling is explained in Appendix C.) Execution can con- 
tinue, however, if the extended error handling facility 
was selected at system generation (FORTRAN 1v (os) 
only). This facility is explained in detail in the 
FORTRAN IV (G and H) Programmer's Guide listed 
in the Preface. 

If a sense light is to be tested and its status recorded, 
entry name SLITET is used. Regardless of its status 
before the test, after a sense light is tested, it is always 
set off. The source language statement is: 


CALL SLITET (i,/) 
14 


where: 
i has a value of 1, 2, 3, or 4, and indicates which 
sense light to test. 
j has a value returned by the subprogram. 1 indic- 
cates the sense light was on; 2 indicates the sense 
light was off. 
If the value of i is not 1, 2, 3, or 4, error message 


216 is issued and execution of this modyle or phase is 


terminated unless the extended error handling facility 


is in effect. 


Overflow Indicator Subprogram 

Entry Name: OVERFL ' 
This subprogram tests for an exponent overflow or 
underflow exception and returns a value that indicates 
the existing condition. After testing, the overflow 
indicator is turned off. This subprogram is called by 
using the entry name OVERFL in a CALL statement. The 
source language statement is: 


CALL OVERFL (/) 


The value of j is returned by the subprogram to indi- 
cate the following: 
1 = floating-point overflow condition has occurred 
Jast. 
2 = no overflow or underflow condition has occurred. 
3 = a floating-point underflow condition has occurred 
last. 


Note: A value for; of 1 or 3 indicates that that condition 
was the last one to occur. An overflow followed by an 
underflow in the same statement would be recorded 
as condition 3 — “underflow occurred last.” 

A detailed description of each exception is given in 
the programmer’s guides listed in the Preface. 


Divide Check Subprogram 
Entry Name: DVCHK 
This subprogram tests for a divide-check exception 
and returns a value that indicates the existing condi- 
tion. After testing, the divide-check indicator is turned 
off. This subprogram is called by using entry name 
DVCHK in a CALL statement. The source language state- 
ment is: 

CALL DVCHK (j) 


where: 
j is set to 1 if the divide-check indicator was on; or 
to 2 if the indicator was off. 


Utility Subprograms 

The utility subprograms perform two operations for 
the FORTRAN programmer;r: they either terminate execu- 
tion (exir) or dump a specified area of storage 
(puMP/PDUMP ). 


End Execution Subprogram 
Entry Name: EXIT 


The end execution subprogram terminates execu- 
tion of the load module or phase and returns control to 
the operating system. (Except that no operator message 
is produced, Exir performs a function similar to that 
performed by the stop statement.) This subprogram 
is called by using the entry name Extr in a CALL state- 
ment. The source language statement is: 


CALL EXIT 


Storage Dump Subprogram 
Entry Names: DUMP/PDUMP 


This subprogram dumps a specified area of storage. 


Either of two entry names (pump or ppuMP) can be used 
to call the subprogram. The entry name is followed by 
the limits of the area to be dumped and the format 
specification. The entry name used in the CALL state- 
ment depends upon the nature of the dump to be 
taken. 

If execution of the load module or phase is. to be 
terminated after the dump is taken, entry name DuMP 
is used. The source language statement is: 
where: 


CALL DUMP (a, b, hy, eves aa; ba, fa) 


@ and b are variables that indicate the limits of stor- 
age to be dumped (either a or b may represent the 
upper or lower limits of storage). 

f indicates the dump format and may be one of the 
integers given in Table 11. The formats available 
depend upon the compiler in use. A sample printout 
for each format is given in Appendix E. 


Table 11, DUMP/PDUMP Format Specifications 


Basic FORTRAN IV FORTRAN IV 


0 specifies hexadecimal 0 specifies hexadecimal 
4 specifies INTEGER 1 specifies LocicaL *1 
5 specifies REAL 2 specifies LOGICAL *4 - 
6 specifies DOUBLE PRECISION | 3 specifies INTEGER *2 

4 specifies INTEGER *4 

5 specifies REAL *4 

6 specifies REAL *8 

7 specifies COMPLEX *8 
8 specifies comPLEx °16 
9 specifies literal 





















If execution is to be resumed after the dump is taken, 
entry name ppuMp is used. The source language state- 
ment is: 


CALL PDUMP (a, b, hy eee, Gny bo, fn) 


where a, b, and f have the same meaning as for DuMp. 


Programming Considerations 

A load module or phase may occupy a different area of 
sturage each time it is executed. To ensure that the 
appropriate areas of storage are dumped, the following 
conventions should be observed. 

If an array and a variable are to be dumped at the 
same time, a separate set of arguments should be used 
for the array and for the variable. The specification 
of limits for the array should be from the first element 
in the array to the last element. For example, assume 
that A is a variable in COMMON, B is a REAL number, 
and TABLE is an array of 20 elements. The following 
call to the storage dump subprogram could be used 
to dump TABLE and B in the hexadecimal format and 
terminate execution after the dump is taken: 


CALL PDUMP (a: bh, hi, o+yQny ba, fo) 


If an area of storage in comMMoON is to be dumped 
at the same time as an area of storage not in COMMON, 
the arguments for the area in comMoN should be given 
separately. For example, the following call to the 
storage dump subprogram could be used to dump the 
variables a and B in REAL*8s format without terminat- 
ing execution: 


CALL PDUMP (A,A,6,B,B,6) 


If variables not in ComMMON are to be dumped, each 
variable must be listed separately in the argument 
list. For example, if k, Pp, and Q are defined implicitly 
in the program, the statement 


CALL PDUMP (R,R,5,P,P,5,0,Q,5) 


should be used to dump the three variables. If the 
statement — 


CALL PDUMP (R,Q,5) 


is used, all main storage between R and Q is dumped, 
which may or may not include Pp, and may include 
other variables. 

If an array and a variable are passed to a subroutine 
as arguments, the arguments in the call to the storage 
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dump subprogram in the subroutine should specify the 
parameters used in the definition of the subroutine. 
For example. if the subroutine sust is defined as: 


SUBROUTINE SUBI (X.Y) 
DIMENSION X (10) 


and the call to sust within the source module is: 


DIMENSION A (10) 


e 
e 


CALL SUBI (A, B) 
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then the following statement should be used in susr to 
dump the variables in hexadecimal] format without ter- 
minating execution: 
CALL PDUMP (X(1), X(J0), 0, Y, Y, 0) 
If the statement 
CALL PDUMP (X(1), Y, 0) 

is used, all storage between a(i1) and y is dumped 
because of the method of transmitting arguments. 

When hexadecimal (0). or literal (9) is specified, 


the programmer should realize that the upper limit 
is assumed to be of length 4. 


Algorithms 


This chapter contains information about the method by which each mathematical 
function is computed. The information for explicitly called subprograms is arranged 
alphabetically according to the specific function of each subprogram (i.e., absolute 
value, exponentiation, logarithmic, etc.). The individual entry names associated 
with each subprogram are arranged logically from simple to complex within each 
function, For example, the heading “Square Root Subprograms” will have algo- 
rithms arranged in the following order by entry name: sQRrT, DSQRT, CSQRT, CDSQRT. 


Information for the implicity called subprograms is arranged alphabetically 
according to function, and alphabetically by entry name within that function. For 
example, the heading “Complex Multiply and Divide Subprograms” will have 
algorithms arranged in the following order: copvp# /cpMpy#, covo#/cMPY#. 


The information for each subprogram is divided into two parts. The first part 
describes the algorithm used; the second part describes the effect of an argument 
error upon the accuracy of the answer returned. 


The presentation of each algorithm is divided into its major computational 
steps; the formulas necessary for each step are supplied. For the sake of brevity, 
the needed constants are normally given only symbolically. (The actual values 
can be found in the assembly listing of the subprograms.) Some of the formulas 
are widely known; those that are not so widely known are derived from more 
common formulas. The process leading from the common formula to the computa- 
tional formula is sketched in enough detail so that the derivation can be recon- 
structed by anyone who has an understanding of college mathematics and access 
to the common texts on numerical analysis! Many approximations were derived 
by the so-called “minimax” methods. The approximation sought by these methods 
can be characterized as follows. Given a function f(x), an interval I, the form of 
the approximation (such as the rational form with specified degrees), and the 
type of error to be minimized (such as the relative error), there is normally a 
unique approximation to f(x) whose maximum error over I is the smallest among 
all possible approximations of the given form. Details of the theory and the 
various methods of deriving such approximation are provided in the reference.’ 
The accuracy figures cited in the algorithm sections are theoretical, and they do 
not take round-off errors into account. Minor programming techniques used to 
minimize round-off errors are not necessarily described here. 


The accuracy of an answer produced by these algorithms is influenced by two 
factors: the performance of the subprogram (see the chapter, “Performance 
Statistics”) and the accuracy of the argument. The effect of an argument error 
upon the accuracy of an answer depends solely upon the mathematical function 
involved and not upon the particular coding used in the subprogram. 


A guide to the propagational effect of argument errors is provided because 
argument errors always influence the accuracy of answers whether the errors are 
accumulated prior to use of the subprogram or introduced by newly converted 
data. This guide (expressed as a simple formula where possible) is intended to 
assist users in assessing the effect of an argument error. 


1 Any of modem numerical anal texts may be used as a reference. One such text is A. Ralston’s A First 
alg in Numerical Analysis (McGraw-Hill Book Company, Inc., New York, 1965). Background information 


algorithms that use continued fractions may be in H. S. Wall’s Analytic Th Continued 
Peter Van Nostrand Co. Inc., Princeton, » Ju. 1948). aire cereale: peat eanene 
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The following symbols are used in this chapter to describe the effect of an 
argument error upon the accuracy of the answer: 


SYMBOL EXPLANATION 
g(x) The result given by the subprogram. 
f(x) The correct result. 
f(x) — g(x) The relative error of the result 
: f(x) given by the subprogram. 
8 The relative error of the argument. 
E |f(x) —e(x)| The absolute error of the result 
| given by the subprogram. 
A The absolute error of the argument. 


The notation used for the continued fractions complies with the specifications 
set by the National Bureau of Standards.? 


Although it is not specifically stated below for each subroutine, the algorithms 
in this chapter were programmed to conform to the following standards governing 
floating-point overflow/ underflow. 

1, Intermediate underflow and overflows are not permitted to occur. This pre- 
vents the printing of irrelevant messages. 

2. Those arguments for which the answer can overflow are excluded from the 
permitted range of the subroutine. This rule does not apply to cps and CABS. - 

3. When the magnitude of the answer is less than 16~®, zero is given as the 
answer. If the floating-point underflow exception mask is on at the time, the 
underflow message will be printed. 


Control of Program Exceptions in Mathematical Functions 


The FORTRAN mathematical functions have been coded with careful control of error 
situations. A result is provided whenever the answer is within the range repre- 
sentable in the floating-point form. In order to be consistent with FORTRAN control 
of exponent overflow/underflow exceptions, the following types of conditions are 
recognized and handled separately. 


When the magnitude of the function value is too large to be represented in the 
floating-point form, the condition is called a terminal overfow; when the magni- 
tude is too small to be represented, a terminal underflow. On the other hand, if the 
function value is representable, but if execution of the chosen algorithm causes an 


overflow or underflow in the process, this condition is called an intermediate over- 
flow or underflow. 


Function subroutines in the ForTRAN library have been coded to observe the 

following rules for these conditions: 

1. Algorithms which can cause an intermediate overflow have been avoided. 
Therefore an intermediate overflow should not occur during the execution 
of a function subroutine of the library. 

2. Intermediate underflows are detected and not allowed to cause an interrupt. 
In other words, spurious underflow signals are not allowed to be given. Com- 
putation of the function value is successfully carried out. 

3. Terminal overflow conditions are screened out by the subroutine. The argu- 


ment is considered out of range for computation and an error diagnostic 
{s given. 


*For'more information, see Milton Abramowitz and Irene A. Stegun (editors), Handbook of Mathematical 
Funttions, Applied Mathematics Serics-55 (National Bureau of Standards, Washington, D.C., 1965). 
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4. Terminal underflow conditions are handled by forcing a floating-point under- 
flow exception. This provides for the detection of underflow in the same 
manner as for an arithmetic statement. Terminal underflows can occur in the 
following function subroutines: EXP, DEXP, ATAN2, DATAN2, ERFC, and DERFC. 


For implicit arithmetic subroutines, these rules do not apply. In this case, both 
terminal overflows and terminal underflows will cause respective floating-point 
exceptions. In addition, in case of complex arithmetic (implicit multiply and 
divide), premature overflow/ underflow is possible when the result of arithmetic 
is very close to an overflow or underflow condition. 
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Explicitly Called Subprograms 


Absolute Value Subprograms 


CABS/CDABS 

1. Write x + iy, = a + ib. 

2. Let vu; = max ( : axl, iy’ ), and bv. = min ( |, | ) 

3. If characteristics of v, and v. differ by 7 (15 for Gas) or more, or if v, = 0, 
thena = v,, b = 0. 

4, Otherwise, 


a=2*u,° yu + %4(=)',andb = 0 
1 


If the answer is greater than 16", the floating-point overflow interruption will 
take place (see Appendix C). The algorithms for both complex absolute value 
subprograms are identical. Each subprogram uses the appropriate real square root 
subprogram (SQRT or DSQRT ). 


Arcsine and Arccosine Subprograms 


ARSIN/ARCOS 
Algorithm 
1. 1f0 SxS ¥, then compute arcsin (x) by a continued fraction of the form: 
arcsin (x) =x + x8° F where 
dt dy 
O° Fay Fa): 
The coefficients of this formula were derived by transforming the minimax 


rational approximation (in relative error, over the range 0x? = %) for 
arcsin (x)/x of the following form: 


arcsin(x) GQ, + dex" 
= = do + x? 3 1 
by + by + bx? + x' +x 
Minimax was taken under the constraint that ay = 1 exactly. The relative error 


of this approximation is less than 2~7*-3, 
IfO =x SS 4, arccos(x) is computed as: 





arccos(x) =-——-— arcsin(rx). 
2 


2. If % < x1, then compute arccos(x) essentially as: 


arccos(x) = 2° arcsin (+ s *), 


This case is now reduced to the first case because within these limits, 


a. jl~*,— 
Os 5 = 








This computation uscs the rea] square root subprogram (sQrr) 
If % <x 1, arcsin(x) is computed as: 


arcsin(x) = ——— arccos(x). 


2 


Implementation of the above algorithms (steps 1 and 2) were carried out with 
care to minimize the round-off errors. 
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3. If — 1x < 0, then arcsin(x) = — arcsin [x| 
and arecos(x) = 2 — arccos |x|. 


This reduces these cases to one of the two positive cases. 


Effect of an Argument Error 


A 


E~ ia For small values of x, E ~ A. Toward the limits (+ 1) of the 


range, a small A causes a substantial error in the answer. For the arcsine, « ~ 6 
if the value of x is small. 





DARSIN/DARCOS 


Algorithm 
1. If0 =x ¥, then compute arcsin(x) by a continued fraction of the form: 
arcsin(x) == x + x3° F where 
d, d, ds d, 
ES ae ee ee 

(x? + co) + (x? +3) + (x? +04) + (x? + Cs) 
The relative error of this approximation is less than 2—57-?, 
The coefficients of this formula were derived by transfo-ming the minimax 


rational approximation (in relative error, over the range 0 x°= %) for. 
aresin(x)/x of the following form: 
arcsin(x) i Q, + aox? + ayxt + ayx® + = | 
x hoe Dye + bot + bx? + 2 | 
Minimax was taken under the constraint that ay) = 1 exactly. 
If 0=x SS ¥, arccos(x) is computed as: 
arccos(x) = — arcesin(x). 


2. 1f %& < x1, then compute arccos (x) essentially as: 


arccos(x) = 2° arcsin (J) ; s), 


This case is now reduced to the first case because within these limits, 


o< Jes». 


This computation uses the real square root subprogram (DsQrT). 
If % < x1, arcsin(x) is computed as: 








arcsin(x) = a — arccos(x). 


Implementation of the above algorithms (steps 1 and 2) were carried out with 
care to minimize the round-off errors. 


3. If - 1 Sx < 0, then arcsin(x) = — arcsin |x|, and arccos(x) = x — arccos |z!. 
This reduces these cases to one of the two positive cases. _ 


Effect of an Argument Error 


A 
E~ Vio * For small values of x, E ~ A. Toward the limits (+ 1) of the 


range a small A causes a substantial error in the answer. For the arcsine, « ~ 6 if 
the value of x is small. 
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Arctangent Subprograms 


ATAN — Basic FORTRAN IV (OS) 


Algorithm 
1. Reduce the computation of arctan (x) to the case 0 S x S 1, by using 


arctan (— x) = — arctan (x), or 


e' 1 — ul D 
arctan { 77 J =-5-— arctan |x|. 


2. If necessary, reduce the computation further to the case |x| < tan 15° by using 


: V3ex-—1 
arctan (x) = 30° + arctan Coy va J 
V3 °“x—- 
The value of os cA <= tan 15° if the value of x is within the range, 
tan 15° <xSl. The value of ( \/ 3° x — 1) is computed as 
(V3 —1)x —1+ x to avoid the loss of significant digits. 
3. For |x| S tan 15°, use the approximation formula: 


rctan 
arctan (2) _ 960310579 — 0.05160454x? 4+ —-0:9913709 


x x* + 1,4087812 ° 
This formula has a relative crror less than 2-7-1 and can be obtained by 
transforming the continued fraction 
x2 


arctan (x) x? 5 


= 


x 3 5 -) 
(2+: — Ww 





where tw has an approximate value of ( - Br + = 10-*, but the true 
4°5 
7°79 


value of w is 


(8 . .\). 
71t* )+ 


The original continued fraction can be obtained by transforming the Taylor 
series into continued fraction form. 


Effect of an Argument Error 


E~ eee For small values of x, « ~ 8; as the value of x increases, the effect 


of § upon e diminishes. 


ATAN/ATAN2 — FORTRAN IV (OS) 


Algorithm 

1. For arctan (2;, x2): 
If x, < 0, use the identity arctan (x;, x2) = —arctan (—x, x2). 
Hence we may assume that x; = 0. Then: 


x 
If either x. = 0 or KE > 274, the answer = c% 
2 


If x2 << Oand Ke < 2-4, the answer = ». 


= ). and 


For the general case, if +. > 0, the answer = arctan (= 
a Xy 
if x.'< 0, the answer = x — arctan aa 


Xe 
x : 
on above, or of arctan(x) for the single argu- 











2. The computation of arctan ( 


ment case, follows the algorithm given for the subprogram aTAN in Basic 
FORTRAN Iv (Os). 


Effect of an Argument Error 


A : 
Ew~ lee For small values of x, « ~ 8; as the value of x increases, the effect 


of § upon ¢ diminishes. 


DATAN — Basic FORTRAN IV (OS) 


Algorithm 
1. Reduce the computation of arctan (x) to the case 0 = x = I by using 
arctan( —~x) = — arctan(x) and 
1 Tv ‘ | 
arctan . , =-——— arctan jx). 
ix] 2, | 


2. If necessary, reduce the computation further to the case |x| S tan 15° by using 


. 3ex—1 
arctan(x) = 30° + arctan (=) 


_ x+ V3 
V3ex-—1 
The ve of aa t ae 
15° < x1. The value of (V/3* x — 1) is computed as (V3 —1)x-14+2% 
to avoid the loss of significant digits. 
3. For |x| < tan 15°, use a continued fraction of the form: 
arctan(x) a; az a3 
cacao <Page 2 ee i a 
x a1 + x| by (by + x7) — (bo + x?) - Ta | 
The relative error of this approximation is less than 2-87, 
The coefficients of this formula were derived by transforming a minimax 
rational approximation (in relative error, over the range 0 S x? S 0.071797) 
for arctan(x)/x of the following form: 
arctan(x) i [3 +,x? + Coxt + Cyx® 
cs ‘3 do + dix? + doxt + x® , 
Minimax was taken under the constraint that ap = 1 exactly. 


SS tan 15°, if the value of x is within the range tan 





Effect of an Argument Error 


Ew~ Ss For small values of x, « ~ 8, and as the value of x increases, the effect 


of « upon § diminishes. 


DATAN/DATAN2 — FORTRAN IV (OS) 


Algorithm — 

1. For arctan(x;, x2): 
If x, < 0, use the identity arctan(x,, x.) = —arctan ( —%;, 2). 
Hence we may assume that x, = 0. Then: 


x 
If either x2 = 0 or He > 258 the answer = oo 
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If x. < O and ey < 2-54 the answer = 7x. 
2 


. 
For the general case, if x, > 0, the answer = arctan (4). and 


. X1 
if x. < 0, the answer = 2 — arctan (3). 


x s 
. The computation of arctan (4) above, or of arctan(x) for the single argu- 
2 


ment case, follows the algorithm given for the subprogram paTAN in Basic 
FORTRAN Iv (0S). 


Effect of an Argument Error 


E~ 


aes For small values of x, ¢ ~ 8, and as the value of x increases, the effect 


of « upon § diminishes. 


Error Functions Subprograms 


ERF/ERFC 
Algorithm 


i 
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If 0 = x S 1, then compute the error function by the following approximation: 
. erf(x) = x(@ 9 + a,x? + aaxt +... + agx!®). 


The coefficients were obtained by the minimax approximation (in relative 
error) of erf(x)/x as a function of x? over the range-0 S x? S 1. The relative 
error of this approximation is less than 2—74-6. The value of the complemented 
error function is computed as erfe(x) = 1 — erf(x). 


. If 1 < x < 2.040452, then compute the complemented error function by the 


following approximation: 

erfc (x) = bo + Byz + bez? +... + boz® 
where z = x — Ty and Ty = 1.709472. The coefficients were obtained by the 
minimax approximation (in absolute error) of the function f(z) = erfe(z + 
Ty.) over the range — 0.709472 = z = 0.33098. The absolute error of this 
approximation is less than 2-°1-5, The limits of this range and the value of the 
origin Ty) were chosen to minimize the hexadecimal round-off errors. The value 


l | 
of the complemented error function within this range is between 356 and 0.1573. 


The value of the error function is computed as erf(x) = 1 — erfc(x). 


. If 2.040452 = x < 13.306, then compute the complemented error function by 


the following approximation: 

erfc(x) = e-*° F/x where z = x? and 
Cy + Co% + C32? 
dyz + doz? + 23° 
The coefficients for F were obtained by transforming a minimax rational 
approximation (in absolute errors, over the range 13.306~? < w S 2.040452-?) 
of the function f(w) = erfc(x) *x* e”, w = x~2, of the following form: 

f(w) _# + ayw + a,w? + ayw* 

oa bo + biw + w? 


The absolute error of this approximation is less than 2~“¢-!, This computation 
uses the real exponential subprogram (Exp). 


F =c) + 


If 2.040452 < x < 3.919206, then the error function is computed as 
erf(x) = 1 — erfe(x). 
If 3.919206 < x, then the error function is = 1. 

4, If 13.306 S x, then the error function is = 1, and the complemented error func- 
tion is = 0 (underflow). 

5. If x < 0, then reduce to a case involving a positive argument by the use of the 
following formulas: 


erf(— x) = — erf(x), and erfe (— x) = 2 — erfe(x). 


Effect of an Argument Error 

.E ~ e+, For the error function, as the magnitude of the argument exceeds 1, 
the effect of an argument error upom the final accuracy diminishes rapidly. For 
small values of x, « ~ 8. For the complemented error function, if the value of x is 


greater than 1, erfc(x) ~ —— Therefore, « ~ 2 x’* 8. If the value of x is negative 


or less than 1, then. ~ e7?’° A, 


DERF/DERFC 

Algorithm 

1. If 0 Sx < 1, then compute the error function by the following approximation: 
. erf(x) S= x(@o + a,x? + aoxt +... + ay1x"?), 


The coefficients were obtained by the minimax approximation (in relative 
error) of erf(x)/x as a function of x° over the range 0 S x? <1. The relative 
error of this approximation is less than 2—5¢-°, The value of the complemented 
error function is computed as erfe(x) = 1 — erf(x). 

2. If 1S x < 2.040452, then compute the complemented error function by the 
following approximation: 


erfe(x) = bo + bz + bez? +... dig} 


where z = x — Ty and Ty = 1.709472. The coefficients were obtained by the 
minimax approximation (in absolute error) of the function f(z) = erfe(z + To) 
over the range —0.709472 <= z <= 0.33098. The absolute error of this approxi- 
mation is less than 2~®-3, The limits of this range and the value of the origin 
To were chosen to minimize the hexadecimal round-off errors. The value of the 


complemented error function within this range is between = and 0.1573. The 


- value of the error function is computed as erf(x) = 1 — erfc(x). 
3.. If 2.040452 = x < 13.306, then compute the complemented error function by 
the following approximation: 


erfc(x) = e~** F/x where z = x” and 
Pi ene... ne, A 
(Z+c1) + (Z+e2) + °°" (Z +06) + (ZH E7)' 

The coefficients for F were derived by transforming a minimax rational approxi- 
~ mation (in absolute errors, over the range 13.306-? = w < 2.040452-?) of the 

function f(w) = erfe(x) + x+e“, w = x~—?, of the following form: 

Qy + ayw + aw? +... + aw? 
(Ola bie tba eo + hae 
The absolute error of this approximation is less than 2~5?-°, This computation 
uses the real exponential subprogram (pexp). If 2.040452 = x < 6. Ogee, then 


the error function is computed as erf(x) = 1 — erfe(x). 
If 6.092368 = x, then the error function is = 1. 


F=Co 
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4. If 13,306 =x, then the error function is = 1, and the complemented error 
function = 0 (underflow). 

5. If x < 0, then reduce to a case involving a positive argument by the use of the 
following formulas: 


erf(— x) = — erf(x), and erfc (— x) = 2 — erfe(x). 


Effect of an Argument Error 


E ~ e~??° A. For the error function, as the magnitude of the argument exceeds 
1, the effect of an argument error upon the final accuracy diminishes rapidly. For 


small values of x, « ~ 8 For the complemented error function, if the value of 
e-*? 


x is greater than 1, erfc(x) ~ on Therefore, « ~ 2x? ¢ 8, If the value of x 
is negative or less than 1, thene ~ e*7 « A. 





Exponential Subprograms 


EXP 
Algorithm 


1. Ifx < — 180.218, then 0 is given as the answer via floating-point underflow. 
2. Otherwise, divide x by log,2 and write 


~_ = 4g-b-—d 
lo 


where a and b are integers,0 = b= 3 and0Sd< 1. 
3. Compute 2~¢ by the following fractional approximation: 


2d 
aE ae 
0.034657359 d? + d + 9,9545948 — ad? + 87.417497 


This formula can be obtained by transforming the Gaussian continued fraction 





y= 


2-622 ] — 


I+ 2-— 3+ 2— 5+ 2% T+. J 
The maximum relative error of this approximation is 2-?°. 
4, Multiply 2-¢ by 2-°. . 
5. Finally, add the hexadecimal exponent a to the characteristic of the answer. 
Effect of an Argument Error 


« ~ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8° x. 


DEXP 
Algorithm 


1. If x < — 180.2187, then 0 is given as the answer via floating-point underflow. 
2. Divide x by log,2 and write 


x= (40 - 5 - 5) + toga — + 
where a, b, and c are integers, 0S b= 3, 0Sc S15, and the remainder r is 


: 1 
within the range OS r < 76° log.2. This reduction is carried out in an extra 
precision to ensure accuracy. Then e* = 16¢ ° 2% « 2-0/8» @—F, 


3. Compute e~" by using a minimax polynomial approximation of degree 6 over 





the range 0 <r < — ° log.2. In obtaining coefficients of this approximation, 


16 
the minimax of relative errors was taken under the constraint that the constant 
term dy shall be exactly 1. The relative error is less than 2—5°-87, 

4. Multiply e-" by 2-°/16, The 16 values of 2-¢/!6 for 0 <= c S 15 are included in 
the subprogram. Then halve the result b times. 

5. Finally, add the hexdecimal exponent of a to the characteristic of the answer. 


Effect of an Argument Error 
E ~ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8° x. 


CEXP/CDEXP 


Algorithm 

The value of e*+ is computed as e” * cos(y) + i* et * sin(y). The algorithms for 
both complex exponential subprograms are identical. Each subprogram uses the 
appropriate real exponential subprogram (Exp or DExP) and the appropriate real 
sine/cosine subprogram (cos/SIN or DCOS/DsIN ). 


Effect of an Argument Error 


The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If e+” = R+ e!¥, then H = yande(R) ~ A (x). 


Gamma and Log Gamma Subprograms 


GAMMA/ALGAMA 


Algorithm 

1. If0 < x S 2-52, then compute log-gamma as log.I'(x) = — log.(x). 
This computation uses the real logarithm subprogram (ALO). 

2. If 2-252 < x < 8, then compute log-gamma by taking the natural logarithm of 
the value obtained for gamma. The computation of gamma depends upon the 
range into which the argument falls. 


Tix +1 
3. If 2-752 < x < 1, then use ['(x) = —— to reduce to the next case. 


4. If 1 SxS 2, then compute gamma by the minimax rational approximation (in 
absolute error) of the following form: 

% [do + 4,Z% + G22? + agz?] 
bo + biz + bez? + 28 


where z = x — 1.5. The absolute error of this approximation is less than 2~25-° 
5. If2 <x < 8, then use P(x) = (x — 1) T(x — 1) to reduce step by step to the 
preceding case. 
6. If 8 = x, then compute log-gamma by the use of Stirling’s formula: 
log.T'(x) = x(log.(x) — 1) — % log.(x) + % log.(2e) + G(x). 
The modifier term G(x) is computed as 
G(x) = dox—! +d,x-2. 
These coefficients were obtained by a form of minimax approximation minimiz- 
ing the ratio of the absolute error to the value of x. The absolute error is less 
than x * 2-62, Remembering the fact that x < log.I'(x) in this range, the 
contribution of this error to the relative error of the value for log-gamma is less 


I(x) eo + 


Algorithms 27 


than 2-282, This computation uses the real logarithm subprogram (ALoc). 
For gamma, compute I'(x) = e”, where y is the value obtained for log-gamma. 
This computation uses the real exponential subprogram (ExP ). 


Effect of an Argument Error 


«e ~ wy (x) ° A for gamma, ae E ~ y (x) ¢ A for log-gamma, where y is the 
digamma function. 


J 
If a << 3, then —2 < y (x) < 1. Therefore, E ~ A for log-gamma. How- 


ever, because x = 1 and x = 2 are zeros of the log-gammna function, even a small 
8 can cause a substantial « in this range. 

If the value of x is large, then y (x) ~ log, (x). Therefore, for gamma, 
« ~ & x log, (x). In this case, even the round-off error of the argument con- 


tributes greatly to the relative error of the answer. For log-gamma with large 
values of x,« ~ 8. 


DGAMMA/DLGAMA 

Algorithm 

1. If0 <x 2-752, then compute log-gamma as log,T'(x) = — log.(x). 
This computation uses the real logarithm subprogram (pDLOc). 

2. If 2-52 < x < 8, then compute log-gamma by taking the natural Jogarithm 
of the value obtained for gamma. The computation of gamma depends upon the 
range into which the argument falls. 


r(ix+1 
3. If 2-*52 < x < 1, then use T(x) = ere to reduce to the next case. 


4. If 1 =x S 2, then compute gamma by the minimax rational approximation (in 
absolute error ) of the following form: 
Z [do + a,z +... + agz*] 
bo + biz +... + bez + zt 
where z = x — 1.5. The absolute error of this approximation is less than 2—5*-8, 
5. If 2 < x < 8, then use r(x) = (x — 1) I(x — 1) to reduce to the preceding 
case. 
6. If 8 = x, then compute log-gamma by the use of Stirling’s formula: 
log.I'(x)} s< x(log.(x) - 1) — % log.(x) + % log.(2r) + G(x). 
The modifier term G(x) is computed as 
G(x) sf dox—! + dyx—3 + dox—5 + dgx-7 + dx, 


These coefficients were obtained by a form of minimax approximation minimiz- 
ing the ratio of the absolute error to the value of x. The absolute error is less 
than x * 2-51, Remembering the fact that x < log,I'(x) in this range, the 
contribution of this error to the relative error of the value for log-gamma is less 
than 2~-5¢!, This computation uses the real logarithm subprogram (pLoc). For 
gamma, compute ['(x) = e”, where y is the value obtained for log-gamma. 
This computation uses the real exponential subprogram (vEXxP). 


T(x) = eo + 


Effect of an Argument Error 


e ~ y(x) ° A for gamma, and E ~ y(x) ° A for log-gamma, where y is the 
digamma function. 


1 
If a <it< 3, then —2 < y(x) < 1. Therefore, E ~ A for log-gamma. How- 


ever, because x = 1 and x = 2 are zeros of the log-gamma function, even a small 
8 can cause a substantial ¢ in this range. 
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If the value of x is large, then y(x) ~ log.{x). Therefore, for gamma, 
« ~ &* x log.(x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of x,<¢ ~ 8. 


Hyperbolic Sine and Cosine Subprograms 


SINH/COSH 


Algorithm 
1. If |xj < 1.0, then compute sinh(x) as: 
sinh(x) Sx + Cyx* + Cox® + yx? 
The coefficient c; were obtained by the minimax approximation (in relative 


ul 


as the function of x2. The maximum relative error of this 





error ) PF 


anacdicaliont 16 2725-8, 
2. Ifx = 1.0, then sinh(x) is computed as: 
sinh(x) = ( 1+ 8) [e+ loge v*/et+ lx.t], 


- , so that this expression is theoretically equivalent to 
[e* — e~*]/2. The value of v (and consequently those of log.t and &) was so 
chosen as to satisfy the following conditions: 

a) vis slightly less than %, so that § > 0 and small. 

b) log.v is an exact multiple of 2-'*. 
The condition b) insures that the addition x + log.v is carried out exactly. This 
maneuver was designed to reduce the round-off errors and also to enlarge the 
limits of acceptable arguments. This computation uses the real exponential sub- 
program (ExP). 

. Ifx S — 10, use sinh(x) = — sinh(|x|) to reduce to case 2 above. 

. If cosh (x) is desired, then for all valid values of arguments use the identity: 
cosh(x) = (1 + 8) [e™+!*" + v?2/e7+!8."], Here the notation and the consid- 
eration are identical to case 2 above. This computation uses the real exponential 
subprogram (ExP),. 


Here, 1 + 8 = 


wm OO 


Effect of an Argument Error 
For the hyperbolic sine, E ~ A * cosh (x) ande ~ A* coth(x). 
For the hyperbolic cosine, E ~ A + sinh(x) ande ~ 8¢ tanh(x). 


Specifically, for the cosine, « ~ A over the entire range; for the sine, e ~ 8 for 
small values of x. 


DSINH/DCOSH 
Algorithm 
1. If |x| < 0.881374, then compute sinh(x) as: 
sinh (x) = Cox + €4x? + Cox +... + cgx)4. 
The coefficients c, were obtained by the minimax approximation (in relative 


rror ) of - = 





as the function of x*. Minimax was taken under the constraint 


that cy = 1 cb: The maximum relative error of this approximation is 2—55-7, 
2. If x = 0.881374, then sinh(x) is computed as: 


sinh(x) = (1 + 8) [e*+!8." — 92/e2+8,0], 
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I 
Here, 1 + 8 = By’ °° that this expression is theoretically equivalent to 


[e? — e-*]/2. The value of v (and consequently those of log.o and 8) was so 
chosen as to satisfy the following conditions: 
a) vis slightly less than 4, so that § > 0 and small. 
b) log.v is an exact multiple of 2~ 1°. 
The condition b) insures that the addition x + loge i is carried out exactly. This 
maneuver was designed to reduce the round-off errors and also to enlarge the 
limits of acceptable arguments. This computation uses the real exponential sub- — 
program (DEXP ). 
3. If x S —0.881374, then use sinh(x) = —sinh(|x|) to reduce to case 2 above. 
4. If cosh(x)is desired, then, for all valid arguments use the identity: 
cosh(x) = (1 + 8) [e7+'8." + v?/et+'."], Here the notation and the consid- 
eration are identical to case 2 above. This computation uses the real exponential 
subprogram (DExP). 
Effect of an Argument Error 
For the hyperbolic sine, E ~ A * cosh(x) ande ~ A* coth(x). 
For the hyperbolic cosine, E ~ 4° sinh(x) ande ~ A* tanh(x). 


Specifically, for the cosine, « ~ A over the entire range; for the sine, ¢ ~ 8 for 
the small values of x. 


Hyperbolic Tangent Subprograms 


TANH. 


Algorithm 
1. If |x| 2-72, then tanh(x) = 
2. If 2-"? < |x| = 0.7, use the following fractional ne ee 
tanh(x) 0.814565 
x? + 2.471749 |" 
The coefficients of this approximation were obtained by taking the minimax 
of relative error, over the range x? < 0.49, of approximations of this form under 


the constraint that the first term shall be exactly 1.0. The maximum relative 
error of this approximation is 2—2¢-4, 


xl]~- | o.o0s78e8 + 


2 
(e*7)? + 1° 
The computation for this case uses the real exponential subprogram (Exp). 
4, If x = 9.011, than tanh(x) = 1. 
5. Ifx < —0.7, then use the identity tanh(x) = —tanh(—x). 


3. If 0.7 <x< 9.01], then use the identity tanh(x) = 1 — 


Effect of an Argument Error 


2A 
E ~ (1 — tanh? x) A, ande ~ sinh (2x) For small values of x, « ~ 8, and as the 
value of x increases, the effect of 8 upon e diminishes. 


DTANH 


Algorithm ; 
1. If |x] S 2-7, then tanh(x) = x. 
2. If Q-VW< || < 0.54931, use the following fractional approximation: 
tanh(x) io d,x? d, ds 
7; §O" Seat Stat Oe o 





This approximation was obtained by rewriting a minimax approximation of the 
following form: 
tanh(x) Ag + yx? + ayxt 
& Cy + 2° ——_-— >. 
by + b,x? + box! +x 
Here the minimax of relative error, over the range x? S 0.30174, was taken 


under the constraint that cy shall be exactly 1.0. The maximum relative error of 
the above is 2- ®. 


3. If 0.54931 S x < 20.101, then use the identity tanh(x+) = 1 — eT 


This computation uses the double precision exponential subprogram (DFxP). 
4, If x = 20.101, then tanh(x) = 1. 
5. Ifx S -- 0.54931, then nse the identity tanh(x) = — tanh(— x). 





Effect of an Argument Error 
2A 
ay = 2 pe eee eta 
E ~ (1 — tanh? x) A, ande sinh (2x) 


value of x increases, the-effect of § upon ¢ diminishes. 


For small values of x, « ~ 5. As the 


Logarithmic Subprograms (Common and Natural) 


ALOG/ALOG10 


Algorithm 


1, Write x = 16? » 2-¢* m where p is the exponent, qg is an integer, 0S q & 3, 
and mis within the range, 4 =m < 1. 
2, Define two constants, a and b ‘(where a = base point and 2-° = a), as follows: 


1 
Ife¥esm<-= Vo’ ,thena = %Zandb =]. 


l 
a= ae al ae 0. 


3. Writez = ae 
m+ 





1 
— and 2 < 0.1716. 


1+ 
=, and log.(x) = (4p — q — b) log. 2 + loa(; = =) 


m—- 


1l+z 
5. To obtain | —_— = Nan 
obtain o8(7— *), first compute w = 2z = 05m + 5 Ea 


resented in our system with slightly more significant digits than z itsetE); and 
apply an psu of the following form: 


(T=) =el-+ 55] 
loge { ———— J = w | co + 
These coefficients were obtained by the minimax rational approximation of 


1 l+z ; 
— 5 tou (7—*) over the range z? « (0, 0.02944) under the constraint that co 


shall be exactly 1.0. The maximum relative error of this approximation is less 
than 2-25-38, 


6. If the common logarithm is desired, then logiox = logice * logex. 


. «Then, m=ae 








1+ 
4, Now, x = 26p—a—-de = 


(which is rep- 





ies of an Argument Error 


~ 8. Specifically, if § is the round-off error of the argument, e.g., 8 ~ 6 * 10-8, 
ae E ~ 6¢ 10-°. Therefore, if the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 
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DLOG/DLOGIO 


Algorithm 


1. Write x = 16 * 2-7 + m where p is the exponent, q is an integer, 0 = q S 3, 
and m is within the range 4% =m <1. 
2. Define two iigiogs a and b (where a = base point and 2—® = a), as follows: 


fen thena = ¥%andb = 1. 








V2’ 
— << = = 
logs m < 1, thena = landb = 0. 
; m— 1 ee ee 
3. Write z = ———.. Then, m = a Dag and |z| < 0.1716. 








1+ | 
4, Now, x = Qfr-e-b 5 and logx = (4p — q — b) log,2 + log, 


m-a we 
5. To obtain log. (5 ) first compute tw = 22 = 05m 4 05a (which is repre- 


sented in our system with slightly more significant digits than z itself), and 
apply an approximation of the following form: 


1+2z C3 
loge (4) s=wl Co t cw? { w? + co + : 
ae te + — 


w? + Ce 
These coefficients were obtained by the minimax rational approximation of 


1 l+z 
m loge i) over the range z? « (0, 0.02944) under the constraint that co 


shall be exactly 1.0. The maximum relative error of this approximation is less 
than 2-69-55, 


6. If the common logarithm is desired, then logiox = logice * log.x. 


Effect of an Argument Error 


E ~ 8. Therefore, if the value of the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 


CLOG/CDLOG 


Algorithm 
1. Write log. (x + iy) = a + ib. 


2. Then, a = log. |x + iy| and b = the principal value of arctan (y, x). 
3. loge |x + iy| is computed as follows: 
Let vo, = max (|x|, |y|), and e2 = min (|x|, |y|). 


1 
Let t be the exponent of 0, i.¢., 0) = me 165, —- Sm. 


> 16 
tift=0 
Finally, let t; = { ' i lift > of. 
and s = 164, 


‘ 0; v2 \? 
Then, log, |x + iy| = 4t, + log.(2) + % log, [ (=) + =) i 


Computation of v;/s and v2/s are carried out by manipulation of the charac- 
teristics of v, and ve. In particular, if v2/s <1, it is taken to be 0. The algor- 
ithms for both complex logarithm-subprograms are identical. Each subprogram 
uses the appropriate real natural logarithm subprogram (Loc or pLoc) and the 
appropriate arctangent subprogram (ATAN2 or DATAN?). 
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Effect of an Argument Error 

The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + iy = r+ e* and log. (x + iy) = a+ ib, thenh = b 
and E(a) = 8(r). 


Sine and Cosine Subprograms 


SIN/COS 
Algorithm 
. 4 
1. Define z = — © |x| and separate z into its integer part (q) and its fraction part 
ms 


-q)+(4-r) 


2. If the cosine is desired, add 2 to g. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of q reduces the general case to the computation 
of sin (x) for x = 0 because 





(r). Then z = q +r, and |x| = (4 


cus (+x) = sin (++ r), and 


sin (— x) = sin(# +x). 
3. Let go = q mod 8. 


Then, for go = 0, sin (x) = sin (+ r), 
qo = 1, sin (x) = cos (Fa - ), 
qo = sn (x) = cos a r), 
qo = 3, sin (x) = sin G a-1)), 
qo = 4, sin (x) = — sin (+--+). 
Go = 5, sin (x) = — cos (Fa _ )), 
qo = 6, sin (x) = — cos (+--+), 
qo = 7, sin (x) = — sin (Fa - n). 
These formulas reduce each case to the computation of either sin (+ ° r) 


or COs (+ ° r) where r; is either r or (1 — r) and is within the range, 
0 = ns 1, 


4, If sin (+ ° ris needed, it is computed by a polynomial of the following 
form: 


sin (+ ° rs) " (do + ayr,? + Query? + asr,° ). 


The coefficients were obtained by the interpolation at the roots of the Chebyshev 
polynomial of degree 4. The relative error is less than 2—?*! for the range. 
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5. If cos (+ °T; ) is needed, it is computed by a polynomial of the following 
form: | 


cos (+ * rm 1+ byr,? + bor,* +bsr,°. 


Coefficients were obtained by a variation of the minimax approximation which 
provides a partial rounding for the short precision computation. The absolute 
error of this approximation is less than 2— 74:57, 


Effect of an Argument Etror 


E ~ A. As the value of x increases, 4 increases. Because the function value dimin- 
ishes periodically, no consistent relative error contro] can be maintained outside 


x wv 
the principal range, — oo Sxs+t > 


DSIN/DCOS 


Algorithm 
1. Divide |x| by and separate the quotient (z) into its integer part (q) and 


4 
its fraction part (r). Then, z = {x| « petit + r, where q is an integer and r 
is within the range, 0 Sr < 1. . 

2. If the cosine is desired, add 2 to q. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of q reduces the general case to the computation of 
sin (x) for x 2 0, because 

cos (+ x) = sin (1 +4), and 
sin (— x) = sin ({x| + 7). 


3. Let go = q mod 8. 


Then, for qo = 0, sin (x) = sin 





(7+) 
go = sin (2) = 00s (F(1~1)), 
qo= % sin (2) = c0s (F-r), 
go = 3,sin (x) = sin (Fa = r)), 
qo = 4,sin (x) = — sin (=-r), 
qo 8,sin (x) = — cs (F(1-2)), 
qo = 6,sin (2) = ~ 00s (+r), 
go = 7,sin (x) = — sin (Za-n). 


nN 





These formulas reduce each case to the computation of either sin (+ ° rn) 


or cos (> “n); where f; is either r or (1 — r), and is within the range, 
0Osnsl. 
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4. Finally, either sin(4 ty Jor cos (+ “ry )is computed, using the polynomial 


interpolations of degree 6 in r,? for the sine, and of degree 7 in r,? for the cosine. 
In either case, the interpolation points were the roots of the Chebyshev poly- 
nomial of one higher degree. The maximum relative error of the sine polynomial 
is 2-58 and that of the cosine polynomial is 2-*-’. 

Effect of an Argument Error 


E ~ 4. As the value of the argument increases, A increases. Because the function 
value diminishes periodically, no consistent relative error control can be main- 


tained outside of the principal range, — zs xs +>. 


CSIN/CCOS 
Algorithm 
1. If the sine is desired, then 
sin(x + iy) = sin(x) * cosh(y) + i* cos(x) ¢ sinh(y). 
If the cosine is desired, then . 
cos(x + iy) = cos (x) * cosh(y) — i* sin(x) ° sinh(y). 
2. The value of sinh(x) is computed within the subprogram as follows. 
Assume x = 0 for this, since sinh( — x) = — sinh(x). 
3. If x = 0.346574, then use sinh (x) = % (« - =; : 
4, If0 =x < 0.346574, then compute sinh(x) by use of a polynomial: 
sinh (x) 


x 





& By + asx? + agx'. 


The coefficients were obtained by the minimax approximation (in relative 
error ) of sinh(x)/x over the range 0 S x? < 0.12011 under the constraint that 


Gy shall be exactly 1.0. The relative error of this approximation is less than 
Q-26.18. 


1 
5. The value of cosh(x) is computed as cosh (x) = sinh|x| + <r 
This computation uses the real expoential subprogram (Exp) and the real 
sine/cosine subprogram (s1n/cos). 
Effect of an Argument Error 


To understand the effect of an argument error upon the accuracy of the answer, the 
programmer must understand the effect of an argument in the szN/cos, Exp, and 
SINH/COSH subprograms. 


CDSIN/CDCOS 
Algorithm 
1. If the sine is desired, then 
sin (x-+iy) = sin(x) * cosh(y) + i* cos(x) « sinh(y). 
If the cosine is desired, then 
cos(x + iy) = cos(x) * cosh(y) — i* sin(x) + sinh(y). 
2. The value of sinh(x) is computed within the subprogram as follows. 
Assume x = 0 for this, since sinh( — x) = — sinh(x). 


3. If x = 0.481212, then use sinh(x) = 4% (e - =): 
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4, 1f0 =x < 0.481212, then compute sinh(x) by use of a polynomial: 
sinh(x) 
x 





St dy $F a,x? + dex? + agx® + ayx® + asx". 


The coefficients were obtained by the minimax approximation (in relative 
error) of sinh(x)/x over the range 0 = x* S 0.23156 under the constraint 


that ay shall be exactly 1.0. The relative error of this approximation is less 
than 2754.07, 


1 
5. The value of cosh(x) is computed as cosh (x) = sinh|z| +o 
This computation uses the real exponential subprogram (pExr) and the real 
sine/cosine subprogram (Dsin/DCos). 
Effect of an Argument Error 


To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument error in the DsIn/DCos, 
DEXP; and psINH/DCOSH subprograms. 


Square Root Subprograms 


SQRT 

Algorithm 

1. If x = 0, then the answer is 0. 

2. Write x = 167°-” + m, where 2p — q is the exponent and q equals cither 0 or 1; 
m is the mantissa and is within the range —> 

3. Then, \’x = 16° + 4-4\/m. 

4, For the first approximation of x, compute the following: 
1.288973 

0.8408065 + m/]° 


This approximation attains the minimax relative error for hyperbolic fits of Vx. 
The maximum relative error is 2- 5-748, 


5. Apply the Newton-Raphson iteration 
Yat = vs (u. +— 
twice. The second iteration is performed as 


n= % (n- =)+4, 


with a partial rounding. The maximum relative error of ye is theoretically 
Q- 25.0. 


d 1. 
i6 = sme< 


Yu = 16° bd 4-4 ° (1.681595 ~" 


Effect of an Argument Error 


1 
e~-z 6. 


2 


DSQRT 


Algorithm 
1. If x = 0, then the answer is 0. 
2. Write x = 16°°-4 + m, where 2p — q is the exponent and q equals either 0 or 1; 


1 
m is the mantissa and is within the range Te sm<l. 
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3. Then, Vx = 16?* 4-7 Vm. 
4. For the first approximation of \/x, compute the féllowiigs 
Yo = 16? © 4'-2 + 0.9202 (m + 0.2587). 


The extrema of relative errors of this approximation for q = 0 are 2-37? at 


1 
m = 1, 2-365 at m = 0.2587, and 2-95 at m = 76° . This approximation, rather 


. x 
than the minimax approximation, was chosen so that the quantity a — y3 be- 
3 


low becomes less than 16°-* in magnitude. This arrangement allows us to 
substitute short form comnterparts for some of the long form instructions in the 
final iteration. 

5. Apply the Newton Raphson iteration 


x 
w= (u +7) 


four times to yo, twice in the short form and twice in the long form. The final 
step is performed as 


x 
Ys= Ys te (= -») 


with an appropriate truncation maneuver to obtain a virtual rounding. The 
maximum relative error of the final result is theoretically 2- &-7%, 


Effect of an Argument Error 


: $ 
e~ 5) 
CSQRT/CDSQRT 
Algorithm 
1. Write Vx + iy = a + ib. 


x| + |x + ty ere eee | 
2. Compute the value 2 = 4] : aERrH, V uy + we where k, w, and we 


are defined in 3, or 4, below. In rs case let v; = max (|x|, |y|) and 


v. = min (|z], |y]). 
3. In the special case when cither v2 = 0 or v; » v3, let w, = v2 and te = v; so 
that w, + we is effectively equal to v,. 
Also me = lif, = |x| and 


= 1/V2if v; = yl. 
4. In the general case, compute F = ‘ | Ye+ V4 (= 2). 
If iilot is near thé underflow threshold, then take 
= |x|, we = 0, ¢ 2F, and k = 1/v2. 
If vo, ¢ i: is near the overflow threshold, then take 
= |x|/4, we = v, *F/2, and k = V2. 


. In all other cases, take w, = [x|/2, we = v, ° F,andk = 1. 
5. Ifz = 0, thena = Oand b = 0. 
If z54 0 and x = 0, thena = z, and 


y 


= oo 


95° 
Ifz>40andx < 0, thena = Le and 
b = (signy) ° z. 
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The algorithms for both complex square root subprograms are identical. 
Each subprogram uses the appropriate real square root subprogram (sQRT or 
DSQRT ). | 


Effect of an Argument Error 


The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + ty =r-e"and Vx + iy = R« ei", 


then «(R) ~ + 8(r), and e(H) ~ 8(h). 


Tangent and Cotangent Subprograms 


TAN/COTAN 
Algorithm 


1. Divide |x! by =z and separate the result into integer part (q) and the fraction 


part (r). Then 'x| = = (qtr). 


2. Obtain the reduced argument (tv) as follows: 
if q is even, then w = r 
if q is odd, then w = 1 — r. 


The range of the reduced argument isO Sw 1. 


Then for qo = 0, tan {x| = tan (+ . w) and cot ‘x; = cot (} . wv), 


= +) 
4 Up 


Go = 2, tan jx; = — cot (+ . v) and cot |x| = — tan (+ . w), 


T 


qv = 1, tan |x| = cot (+ . w) and cot jx| = tan ( 





qo = 3, tan |x| = — tan (+ . w) and cot |x| = — cot (+ . w). 


4. The value of tan (+ ° w) and cot (+ ° w) are computed as the ratio of two 
polynomials: 


v w °P(u) w Q(u) 
tan (F-w =O (4-¢) =a S50 
where u = 4%wu* and 


P(u) = — 8.460901 + u 
Q(u) = — 10.772754 + 5.703366 » u — 0.159321 + u2. 


These coefficients were obtained by the minimax rational approximation (in 
relative error) of the indicated form. The maximum relative error of this 
approximation is 2-26, Choice of u rather than w? as the variable for P and Q 
is to improve the round-off quality of the coefficients. 

. Ifx < 0, then tan(x) = — tan |x|, and cot(x) = — cot |x]. 

. This program is provided with two kinds of error controls, One is for arguments 
whose magnitude is greater than 2!® + x, The other is for arguments which are 
very close to a singularity of the function. In either case, the precision of the 
argument is deemed insufficient for obtaining a reliable result. More specifically, 
the second control screens out the following arguments: 

a) |x| S 16—® for coran (the result would overflow ). 
b) x is such that one can find a singularity within eight units of the last digit 


a 
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value of the floating-point representation of the sum q + r. Singularities are 
cases when the cotangent ratio is to be taken and w = 0. . 
The test threshold of this control can be dynamically modified by assembler 
code programs. | 


Effect of an Argument Error 


A 2 ; Ae 
Ew~ Cos*(x) ,ande ~ sin(ax) for tan(x). Therefore, near the singularities 


1 ee 
x= (: ++) a, where k is an integer, no error control can be maintained. This 


is also true for cotan(x) for x near kr, where k is an integer. 


DTAN/DCOTAN: 
Algorithm 


1, Divide |x| by~, and separate the result into integer part (q) and the fraction 


part (r). Then |x| = (4 +r), 


2. Obtain the reduced argument (tv) as follows: 
if q is even, then wv = r 
if q is odd, then w = 1 — r. 
The range of the reduced argument is 0S w = 1. 
3. Let go = q mod 4. 


Then for qo = 0, tan |x| = tan (+ . w and cot {x| = cot (+ . w), 
Qo = 1, tan |x] = cot (+ ° w )and cot .xj = tan (+ ° wv), 
qo = 2, tan |x| = — cot (+ ° w and cot |x| = — tan (+ . wv), 
Go = 3, tan |x| = — tan (+ . w )and cot ;z| = — cot (+ ° wv). 


4, The value of tan (+ . w) and cot (+ . w) are computed as the ratio of 


two polynomials: 


. (5 ) we P(w?) F (5 ) O(w?) 
n 4 a Nas O(w*) » and cot ra wu = tp P( we)’ 
where both P and Q are polynomials of degree 3 in w. The coefficients of P 
and Q were obtained by the minimax rational approximation (in relative error) 


1 
of mm tan (+ w ) of the indicated form. The maximum relative error of this 


approximation is 2—55-°, 

5. If x < 0, then tan(x) = — tan |x|, and cot(x) = — cot |z|. 

6. This program is provided with two kinds of error controls. One is for argu- 
ments whose magnitude is greater than 2°° « x. The other is for arguments which 
are very close to a singularity of the function. In either case, the precision of 
the argument is deemed insufficient for obtaining a reliable result. More 
specifically, the second control screens out the following arguments: 

a) |x| = 16-° for coran (the result would 0. flow). 

b) x is such that one can find a singularity v :in eight units of the last digit 
value of the floating-point representation o the sum q + r. Singularities are 
cases when the cotangent ratio is to be take and w = 0. 
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The test threslold of this control can be dynamically modified by assembler 
code programs. 


Effect of an Argument Error 


2 
E ~ 


4 
cos*(xy’ and « ~ Sin(Gx) for tan(x). Therefore, near the singularities of 


] . 
x= (: + =) a, where k is an integer, no error control can be maintained. 


This is also true for cotan(x) for values of x near kz, where k is an integer. 


Implicitly Called Subprograms 


The entry point names of the following implicitly called subprograms are gener- 
ated by the compiler. 


Complex Multiply and Divide Subprograms 


CDVD#/CMPY + (Divide/Multiply for COMPLEX*8 Arguments) 
CDDVD#:/CDMPY+ (Divide/Multiply for COMPLEX*16 Arguments) 
Algorithm 
Multiply: (A + Bi) (C + Di) = (AC — BD) + (AD + BC)i 
Divide: (A + Bi)/(C + Di) 
1. If |C| = |D|, set 
A=B,B= —A,C = D,D= —-C,since 
A+ Bi B= Ai 


te ee | DOE tep 2. 
CrD D=Gt 


,_A a B __D 
2. Set A =o:B = Ed HG 
then compute 
A + Bi A’ + Bi A’ + B/D B’— A’D' . 
C+Di~ 1+bDi~ 1+pDbD *t T+pDD"* 


Error Conditions 


Partial underflows can occur in preparing the answer. 


Complex Exponentiation Subprograms | 


FCDXI# (COMPLEX* 16 Arguments) 

FCXPI# (COMPLEX*8 Arguments) 

Algorithm 

The value a Yi + Yyei = (2%, + Z2i)/ is computed as follows. 


Let |j| =) r, ° 2* where r, = U or 1 fork = 0,1,..., K. 
k=0 : 


Then z [4] = Tz, and the factors z** can be obtained by successive squaring. 
, # 


More specifically: 
1, Initially: k = 0, n) = [jf y, + y.4 = 1 + 04, 
a1 (0) Zn (Of = 2% + Zi. 
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9. Raise the index k by 1, and let n'*-") = 2q + 1, where q is the integer 
quotient and r = Oor }. ; 
3. Letn™ = q. 
4. Ifr = 0, then y,® + yo? = yyF—Y + yoo’ 
Ifr = 1, then y,™ + yoi=(yi 8) + yo Yi) (a FOTY Fb ze). 
5. Ifn® +: 0, then!) + 22 = (z,%-D + 34%) )2, and steps 2 
through 5 are repeated until n™ = 0. 
6. When n™) = 0, andj = 0, then y, + yot = yr +y2"i. 
Ifj < 0,theny, + yoi = (1+ 08) / (yi + yo). 


Exponentiation of a Real Base to a Real Power Subprograms 


FDXPD+: (REAL*8 Arguments) 
FRXPR+ (REAL*4 Arguments) 


Algorithm 

1. Ifa = Oand b & 9, error return. 
If a = Oand b > O, the answer is 0. 

2. If a 540 and b = 0, the answer is 1. 

3. All other cases, compute a” as e? + '™ ", In this computation the exponential sub- 
routine and the natural logarithm subroutine are used. If a is negative or if 
b + log a is too large, an error return is given by one of these subroutines. 


Error estimate 

The relative error of the answer can be expressed as (¢. + «) b* log (a). + « 
where ¢;, «2, and ¢; are relative errors of the logarithmic routine, machine multi- 
plication, and the exponential routine, respectively. 

For Fpxpp*, « S 3.5x10~-'8, «. S 2.2x10-'*, and e; SS 2.0x10~-'*, Hence the 
relative error S 5.7x10~1%x | b + log a | + 2.0x10~1°, Note that b + log a is the 
natural logarithm of the answer. 

For Faxpr*, ¢, = 8.3x10-7, «. = 9,.5x10—7, and e3 S 4.7x10-7. Hence the relative 
‘error S 1.8x10~-* x | be loga| + 4.7x10—*. 


Effect of an Argument Error 


[a(1 + 8)] b(1 + &) s= a°(1 + 8,5 loga + b8,). Note that if the answer does 
not overflow, |b * log a| < 175. On the other hand b can be very large without 
causing an overflow of a° if log a is very small. Thus, if a = 1 and if b is very 
large, then the effect of the perturbation 8, of a shows very heavily in the relative 
error of the answer. 


Exponentiation of a Real Base to an Integer Power Subprograms 


FDXPI+ (REAL*8 Arguments) 
FRXPI#t (REAL*4 Arguments) 
Algorithm 


K 
The value of y = a! is computed as follows: Let |j| = )) 1.2* where r;, = Oor 1 


fork = 0,1, ..., K. Then all = 7 a and the factors a** can be obtained by suc- 
r, : 
cessive squaring. 
More specifically: 
1. Initially: k = 0,n = |j{, y = 1, and z = a. 
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2. Raise the index k by 1, and decompose n*~1) = 2q + r, where q is the 
integer quotient and r = 0 or 1. 

3. Let n™ = q. 

4, Ifr = 0, theny™ = y@-0, 
Ifr = 1, then y®) = ye—1ztk—-0), 

5. If n™® ot 0, then z) = z(*-1z(1—) and steps 2 through 5 are repeated 
until n™) = 0, 


1 
6. When n = 0, andj = 0, then y = y). If} < 0, theny = ros 
Note: The negative exponent is computed by taking the reciprocal of the posi- 
tive powey. Thus it is not possible to compute 16.0**—64 because there is a lack 
of symmetry for real floating-point numbers — i.e., 16.0** —64 can be represented, 


but 16.0**64 cannot. The result is obtained by successive multiplications and is 
exact only if the answer contains less than 14 significant hexadecimal digits. 


Exponentiation of an Integer Base to an Integer Power Subprogram 


FIXPI#: (INTEGER*4 Arguments) 
Algorithm 


K, 
The value of L = PF) is computed as follows: Let j = y t, ° 2* where = 9 
k= 0 
or 1 fork = 0,1,..., K. Then I! = z I”, and the factors I** can be obtained by 
t, % 0 
successive squaring. : 
More specifically: 


1. Initially: k = 0,n = j,y = 1, and m® = I, : 

2. Raise the index k by 1, and decompose n‘*~)) = 2q + r, where q is the 
integer quotient and r = Oor 1. 

3. Letn™ = q, | 

4, Ifr = 0, theny™ = y@—3), 
Ifr = 1, theny™ = y—2) 6 muk—1), 

5. If n® at 0, then m'*) = m(t—2) » m(*—)), and steps 2 through 5 are repeated 
until n™ = 0, 

6. When n™ = 0,L = L“), 


Note: The result is obtained by successive multiplications. The result is exact 
only if it is less than (2**31) — 1. Results are meaningless when this limit is 
exceeded and may even be of changed sign. 


This chapter contains accuracy and timing statistics for 
the explicitly called mathematical subprograms. These 
statistics are presented in Tables 12 and 13 and are 
arranged in alphabetical order, according to the entry 
names. The following information is given in the two 
tables: 


Entry Name: This column gives the entry name that 
must be used to call the subprogram. 


Argument Range: This column gives the argument 
range used to obtain the accuracy figures. For each 
function, accuracy figures are given for one or more 
representative segments within the valid argument 
range. In each case, the figures given are the most 
meaningful to the function and range under consid- 
eration. 


The maximum relative error and standard deviation 


of the relative error are generally useful and revealing 
statistics; however, they are useless for the range of a 
function where its value becomes 0, because the slight- 
est error in the argument can cause an unpredictable 
fluctuation in the magnitude of the answer. When a 
small argument error would have this effect, the maxi- 
mum absolute error and standard deviation of the 
absolute error are given for the range. For example, 
absolute error is given for sin(x ) for values of x near =. 


Sample: This column indicates the type of sample 
used for the accuracy figures. The type of sample de- 
pends upon the function and range under consider- 
ation. The statistics may be based either upon an 
exponentially distributed (E) argument sample or a 
uniformly distributed (U) argument sample. 


Accuracy Figures: This column gives accuracy fig- 


ures for one or more representative segments within. 


the valid argument range. The accuracy figures sup- 
plied are based upon the assumption that the argu- 
ments are perfect (i.e., without error and, therefore, 
having no error propagation effect upon the answers). 
The only errors in the answers are those introduced by 
the subprograms, The chapter, “Algorithms,” contains 
a description of some of the symbols used in this 
chapter; the following additional symbols are used in 
the presentation of accuracy figures: 


Performance Statistics 


The maximum 
relative error 
produced 
during testing. 


f(x) — g(x) 


M(e) = Max E(x) 








The standard 
f(x) — g(a) = gin deviation (root- 
F(x) x;) mean-square ) of 
the relative error. 





The maximum 
absolute error . 


produce 
during testing. 


The standard 
deviation (root-- 
a (E) = Ei (x1) - — g(x) yf mean-square ) of 


the absolute error. 


M(E) = Max | f(x) — g(x) | 


In case of complex functions, the absolute value 
signs employed in the above definitions are to mean the 
complex absolute values. In the formulas for the 
standard deviation, N represents the total number of 
arguments in the sample; i is a subscript that varies 
from 1 to N. 


Average Speed: The average time given for each 
function was determined by executing a job using an 
interval timer to measure time differences. Because of 
the various methods of performance enhancement used 
in some models of the System/360 (e.g., interleaving 
storage, buffered storage, execution overlap, antici- 
patory branching analysis), instruction order and mix 
can affect performance on certain models. Conse- 
quently, the times given must be considered as highly 
generalized and differences may be observed during 
normal use of the functions. 

Test ranges, where they do not cover the entire 
legal range of a subroutine, were selected so that users 
may infer from the accuracy figures presented the 
trend of errors as an argument moves away from the 
principal range. The accuracy of the answer deteri- 
orates substantially as the argument approaches the 
limit of the permitted range in several of the sub- 
routines. This is particularly true for trigonometric 
functions. An error generated by any of these subrou- 
tines, however, is at worst comparable in order of 
magnitude to the effect of the inherent rounding error 
of the argument. 
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Table 12. Accuracy Figures 


Accuracy Figures _ 
Entry Argument Sample #+#——---—--—---—____ 
Name Range E/U Relative Absolute 


M (e) o (e) M (E) ¢ (E) 
U 1.16 X 10-6 3.54 x 10°? 


1.25 x 10° 3.04 x 1077 


8.0<X < 160 LiBxlo® 3.80 x 1077 
16.0 =X < 500.0 9.85 X 1077 


| 190x107 | 
055X515 6.85 x 108 
ALOG 10 [055X135 U 7 7.13 x 10° 
1.05 x 10° 2.17 X 107? 
ARCOS 8.85 x 1077 3.19 x 1077 
| ARSIN 9.34 x 10-7 
ATAN The full range 1.01 x 10-* | 
ATAN 2 The full range / 1.01 x 1 


CABS The full range j | Note 1 | 9.15 x 10°? 


Ccos {X;| S 10, |X| 1 250 X 107° - 7.66 X 107 
See Note 2 


U 


CDCOS [Xi; = 10, {X:! S1 U 3.98 x 10°'* 2.50 x 10°tS 
See Note 3 : 





ALGAMA 








ALOG 





2.26 x 10-8 



































1.10 x 107?* 
9.64 x 10°° 


5.38 x 10°}? 
2.25 X 10716 


4.06 x 10717 


2.67 xX 10°? 
2.73 X 10-7 


CDEXP [Xi] S 1, |X| S «/2 
[X:| S 20, |X.| S 20 


The full range 
except (1 + Oi) 


-{Xij S 10, |X| S1 


CDSQRT The full range 


CEXP X,| S 170, |Xel S w/2 
[Xa] = 170, ° U 
«/2 <= |Xa| = 20 
CLOG The full range 
except (1 + Oi) 
-1l0SX <0, U 
ec X10 
emis | vf 
COTAN | [Xi< 4/4 1.07 x 10-8 | 358x107 


V4 < KIS 4/2 
i 


x/2 < |X| S10 1.30 X 10° (Note5) | 311x107 | 
10 < 1s 100 y[sisxi07 | 


Notes: (See end of table.) 














3.76 X 10°! | 
2.74 X 10715 
2.72 x 10716 












CDLOG 





CDSIN 2.35 x 10735 


See Note 4 
1.76 x 107* 
9.93 x 1077 
1.07 x 10° 


7.15 x 10°? 
















4.60 x 10° 
1.28 x 1077 


4.60 x 107° 









1.49 x 10-8 (Note 5 





rs 
mh 


Table 12. Accuracy Figures (Continued) 

















Accuracy Figures 


Oe ee A ore en a a eC 


Relative 
o (e) 


7.38 < 1077 















Sample 
E/U 


Argument 
Range 


Entry 


Name 
1.92 x 10° 


CSIN \X;| S10, |X:| $1 a 
See Note 6 

CSQRT 7.00 x 10" mmxior | 

U_[zorxio | tosxio | 

Ix u_ | 2. ae 

a, 

Ppseeesnes al 





M (e) 











j=l 


7.04 x 10737 


DATAN The full range Note7 | 2.18 x 10°16 
DATAN2 | The full range Note 7 | 2.18 x 107% 7.04 x 10717 


U 
-10=X <0, U | 1.75 x 1076 
U 2.64 x 10735 
9.05 X 10° 


ro XS10 
3.63 X 1071 v7 
2.46 X 10-'* (Note 5 
14 
14 














6.53 x 107? 
5.93 x 10717 









DCOSH | |xX|S5 
2.78 X 1078 (Note 5 


DCOTAN | |X| S x/4 
5.40 X 1073 (Note 5)| 1.13 x 107 


r/4 << |X| Sx/2 
8.61 x 107! (Note5)| 4.61 x 10° 


1.89 X 10°16 2.60 x 1071" 
2.87 X 10717 9.84 x 10738 


~~ | 







1.39 x 10717 8.02 X 10718 
2.08 X 10728 | 652X 1077 | 
| 865x107 


3.26 < 10°%° 8.65 x 10° 


3.51 x 10715 1.96 x 10° 


DGAMMA [0 <X<1 


16=X < 57 


DLGAMA |O0<XS05_ 


2.03 x 10°18 


ré 
7 
7 
: 
7 
7 
7 
8 
7 
6 
7 


wv) 
5 
Q 
fe) 
n 


1 
1 
a 
1 
1 





2.21 x 10736 8.49 x 10° 
5.05 x 10°16 1.90 x 10° 
6.02 x 10715 1.78 X 107'5 
2.77 x 10-31% 9.75 X 10° 
2.89 x 10°36 8.80 x 1071? 
2.86 X 10°'¢ 8.92 xX 10737 


1 
1 
1 
5.43 x 107} 
1 
1 


g 


end of table.) 


4.87 X 10° 
7.77 x 10717 


4.60 x 10°37 2.09 x 107?" 









Norges: (See 
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Table 12, Accuracy Figures (Continued) 


Argument Sample 
Range E 
ange /U M (E) o (E) | 


X<05,X>15 | E | 302x107 6.65 X 1077 
17 


DSIN IX} S «/2 3.60 x 1077 4.82 x 10717 7.74 x 1077 
r/2 < |X] S10 a 1.64 x 10726 
10 < |X| < 100 
IX! < 0.88137 


DSINH 2.06 x 10-76 
all 0.88137 < |X| <5 
The full range 


1.06 x 10°16 
DTAN |X] S x/4 
r/4< |X| S x/2 
r/2 < |X| S10 
10 < |X| S 100 
DTANH | |X| $ 0.54931 
a 0.54931 < |X| $5 


ERF iXj S 1.0 
1.0 < |X| = 2.04 
2.04 < |X| = 3.9192 


-38 << X< A 
0<XS10 















Accuracy Figures 





Absolute 





1.07 X 10727 


1.98 x 10°}? 
6.49 < 10°27 
1.03 x 10°75 










2.16 xX 10737 


6.27 X 10717 
1.43 x 10"? (Note5)} 2.95 x10 | 
2.78 X 10°!3(Note 5) | 7.23 x 10°28 
3.79 x 10°12 (Note 5)| 9.50 X 10-4 
1.91 x 1078 3.86 X 1077 
1.54 x 1076 1.87 x 10-17 
8.16 x 1077 
1.13 x 10°? 
5.95 X 10°8 





















3.70 x 10°® 


9.10 x 10°? 2.96 x 10-7 
1.54 x 10° 3.78 X 10-7 
2.28 x 1076 3.70 x 1077 
1.55 x 1075 8.57 x 1076 


4.65 x 10-7 
9.86 x 10°? 3.66 x 10°? 


1,13 x 1077 
9.47 x 10-7 3.79 Xx 10-7 
2.28 x 107 | 8.32 x 107 
2.20 x 10-8 7.61 x 10-8 







10< Xs 2.04 
204<XK = 40 
40<X5 13.3 
x|S1 
1<|X|/S170 
GAMMA J0<X<10 
10SXS2.0 
20<XS4.0 
40<¢X < 8.0 
80=X S 16.0 


—_— 








16.0 < X $57.0 4.62. X 10-8 1.51 X 10° 

SIN IX| < w/2 1.32 x 10-8 | 118X107 | 
eee ar ee 

. Se iio 

fore [etme | ewww Lownie 


Notes: (See end of table) 
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Table 12. Accuracy Figures (Continued ) 


Entry 
Name 


TAN 













Accuracy Figures 

Absolute 

EC OC LLC) 
7 








BE 2 
ac eee 0 _[A0sxio* (Noes) [asoxie? 
ve leant v[axi?_faasxiee | 
NOTEs: 


* The distribution of sample arguments upon which these statistics are based is exponential radially and is uniform around the 
origin. 
* The maximum relative error cited for the ccos function is based upon a set of 2000 random arguments within the range. In 











the immediate proximity of the points (n + +) x + Oi (where n = 0, + 1, + 2,...,) the relative error can be quite 


2 

high, although the absolute error is small. 
* The maximum relative error cited for the cpcos function is based upon a set of 1500 random arguments within the range. 
In the immediately proximity of the points ( n+ + 


high, although the absolute error is small. 


* The maximum relative error cited for the cpsin function is based upon a set of 1500 random arguments within the range. 
In the immediate proximity of the points nx + Oi (where n = + 1, + 2,...,) the relative error can be quite high, although 
the absolute error is small. 

* The figures cited as the maximum relative errors.are those encountered in a sample of 2500 random arguments within the 
respective ranges. See the appropriate section in the chapter “Algorithms” for a description of the behavior of errors when 
the argument is near a singularity or a zero of the function. A) NE 

* The maximum relative error cited for the cstv function is based upon a set of 2000 random arguments within the range. In _ 

the immediate proximity of the points nr + Oi (wheren = + 1, + 2,...,) the relative error can be quite high, although 


x + Oi (where n = 0, + 1, + 2,...,) the relative error can be quite 











the absolute error is small. 


" The sample arguments were tangents of numbers uniformly distributed between — Sand + =a 


Performance Statistics 47 


-Table 13. Average Machine Timings (Part 1 of 4) 








ATAN2 The full range Note 2 452 
The full range Note 1 116 
pesto st | 0 _| ves 
CDABS The full range Note 1 17009 4628 774 157 
peizia eet |v 08 | 
|X:| S 1, [Xe] S w/2 
|X; S 20, |X:| S 20 


The full range 
(except (0 + Oi) ) 


\X,] = 10, |X! $1 
CDSQRT The full range 
~— TEXy S170, [Xe S 2/2 
[X:] S170, r/2<'X.| < 20 


m= 
be 
iv] 
or 


po] wo 
= 
font 
~] 
a 
an 








Timing Averages in Microseconds for Various Models 
Entry Argument Sample 75 85 

pee | fer Te Pe |e fetal oto | wm 
wigan [OSSE SIO o_| ans | onan [ror] ave] os | a os | w]e) 
[omer ae om fae ie ap sep Da 
aoe% <0 U_ | ioeio [ones | av | —oie| was |e [or | oo | 10, 
See U [tise | avre-| see] ae] oo | moo |e |e 
noes ust me een on me oe 
moe fossxais fw | on nor] wef os] 6 | fe | fT 
al ee os eae 
eee te ee ee eet eat ee 
ion | ee | aos [os [ao [ar | ae | 30 
Toe ae ie [ar [aes so oe 
ATAN TN EN A 
ae = proce | on | a [ es er[ oa [et po |_| 
Note 2 | tos {58 | 2 | | io | 7 
| Note 1 | 39 | 3 
| 35 | 
zz 





os 
v_ [sas 


Note1 | 46248 | 


2 
nS 
ay 
© 





U 60406 | 13995 
Note I 30436 9217 


U 14720 4046 
U 15017 4138 


3596 


1396 


1250 


4 R 
= 3/8 
hr 

— 

oo 

wp 

xs 

a 


rach Nac fe 
| m| & 


— 








The full range 
(except (0 + OF)) 


OSX 


-10=X<0, 
ro XH10 


10<|X}<100 


— | 
fom Be Co) 
Gait 
Got or 


3 
~] 
mF 19 aly 
“1 a1 ul 
re) 
oo 


on 
= 


oo 
NSS 
> 
~] 
oo 


4529 
6399 | 2011 | | 


oo 
> 
vn 





or 
oO 
a 
— 
tt 
a 


f 
—_ rl o Ei 
— 


lals|ale 


ol @ 
als 
o 
~l 
-— 
& 
5 
=] 


IX! S #/4 

w/4 < |X| S /2 
x/2<[X} = 10 
10 < jX| = 100 
IX:} S 10, iX.| = 1 


068 1311 391 
1378 397 


1378 





397 


CSIN 141 


CSQRT 
DARCOS 


NoTEs: (See end of table.) - 


|elelelel|6| « 


EE 





The full range 


Note 1 10446 


U 19127 


Dm | me: 
o 
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Table 13. Average Machine Timings (Part 2 of 4) 


Entry 


Name 


DARSIN 
DATAN 


DATAN2 
DCOS 


DCOSH 
DCOTAN 


DERF 


DERFC 


DGAMMA 


DLGAMA 





S 
8 






Timing Averages in Microseconds for Various Models 
8 
4) 91 195 


Argument Sample [-——— 


75 
Range E/U (Note 3) | (N 


° 
ara 









m1 OO 


The full range 








































20 < |X| = 170 


Ue | 29 | 
og 
10 << X 5100 U | 29 | 17 | 
U 29 | 19° 8 
x/4< |X| S 4/2 U 28 8 
7/2 < |X| S10 U 29 
; S o 
E pee pes 
1.0 < |x| =2.04__ U a7 | 20 | 
2.04 < |X| < 6,092 U 
i ei 
<a rae [es 
a [se [ ase 0) 
v cacweae 
30 caCEE: 
0 ciate 
i<piSm 3609 [a [es [8 
Ce sa fa 


No 
bme 
to 
~— 
i) 
xz 


112 


26290 218} 125 


30931 6706 1469 265 145 
10833 


| 2189 | 
10899 | 2195 | 425| 238 
| 1022 


— 
= 

Gi Ei 

Gi w 
| 


~ 
—~ 


to 
rs 
S| 
Crt 
o> 
co 
a 
— 
~ 
be 
B 
N 
So 
a 
wo] ve 
— | = 
Ga 


> 
G2 
iw) 
BE 


1X < 57 
05<X< 3 
165 X < 500 
05S XS15 
X < 05,X > 15 


05S3XS15 
X < 05,X > 15 


Go 
bo 


GQ} & we. 7 bo 
Go — 


de 


pent 

-_ 
so 
— 





78 
72 


| 


oS 


qj 
a 
S 
2 
~2 
oo 
3 
— 
he 
aw 
BG 
& 


~~ 
he 
hes 
8 
i?) 
a 


Notes: (See end of table. ) 


Performance Statistics 4 


© 


Table 13. Average Machine Timings ( Part 3 of 4) 

















= | et 


i) 
wL 
>} 
“2 


= 8 
=?) 
oD 
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Timing Averages in Meepecen for vers Models 





ria 
ae 
rissa 
| 3598 
ae 
Pia [3 
1124 
1764 
1051 


: 


3 





2679 
ae 
i saeee lise 
eT 
foes |e 
fae 
9020 
9300 
4129 
4633 
4493 
6694 
11729 
11729 
4383 
45 
4560 
3162 
ae 
5378 
|_ S205 | 1364 
a 














nN 
és wn ae an < = es : el. 
o 
Sg elel (Mio| [Fleis|iim) iglele) [SiSi3] lelelelelelSiB] i7/Sielel jFialg 
E¢ =(5/%/S! lywlwitl|3|% wy Vie} St wl Sls jalt jolly Vil FS Vil} Vit ne 
oS 1M ivlsly| & WEST] om et oe] a] w]e | vu vi VOL VE vil se fox facd ML el ele lode ML 
Bm lal llBlV| Els ells] BTV lel] ET Ol lee | [oe] tl Vi elo] Le foe |S] Elz] ef] 2] A! Ele] 
SSS Ey lEISl si siCivd Mist) Vid ol ZL Ly Vwi Yisiwleimilvivig | 3 
V Oe VIV o V ra V Vj Vi vil} Vi VV Vil V Vic 
Vil V {MU eo [MY ede | V bad Dod ie Slsivileo|SlelslViVlelelelelei<i& mL eTSI< Slo] Ss 
eISelRISIEIEISISlslElelEIS Sl Tle lS S/S lz Alol(SlalSl(siS/ELeiS] [Rye ELE Ss] 
as ce ast o 
3 
A 
Pe fe He 3 a 
< be 
ae ae 2 | 3 e |2 52 5 
Ble Ee % < oO} < z 
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Table 13. Average Machine Timings (Part 4 of 4) 


Entry — Arwunent Sample Timing Averages in Microseconds for Various Models 
a cd 


arenes bib ee barber pa ete tet 


1 The distribution of sample arguments upon which these statistics are based is exponential radially and is uniform around the 
origin. 


NOTEs: 


* The sample arguments were tangents of numbers uniformly distributed between — 


* The statistics for the Model 75 are based upon four-way interleaving. 


* The second column of speeds for the Model 85 applies to that machine with high-speed multiply feature. 
* Timing figures on this line apply to the version for Basic FORTRAN IV (os). 


T r 
“7 and + = 
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Appendix A: Assembler Language Information 


The mathematical and service subprograms in’ the 
FORTRAN IW library can be used by the assembler 
language programmer. Successful use depends on three 
things: (1) making the library available to the linkage 
editor; (2) setting up proper calling sequences, based 
on either a call macro instruction or a branch; and (3) 
supplying correct parameters—i.c., arguments. . 


Library Availability 
The System/360 Operating System ForTRAN Iv library 
is a partitioned data sct named sysi.Forriip. The as- 
sembler language programmer must arrange for the 
desired subprograms (modules) to be taken from this 
library and brought into main storage, usually as a part 
of his load module. This can be done by emploving 
the techniques described in the publication IBM 
System/360 Operating System: Linkage Editor and 
Loader, Form C28-6538. : 

For example, the Fortran rv library could be made 
part of the automatic call library by using these job 
control statements: 


:fjobname JOB desired operands 
/Istepname EXEC ASMEPCLG, PARM.LKED 
=’XREF, LIST, 
MAI’ : 
./ ASM. SYSIN Db - 


(assembler language program souree deck) 
7* 
i/LKED SYSLIB DD 
/* 
Subprograms requested in the source program would 


then be available to the linkage editor for inclusion in 
the load module. 


DSNAME=SYS1.FORTLIB, 
DISP=SHR 


Calling Sequences 
Two general methods of calling are possible: (1) cod- 
ing an appropriate macro instruction (see the publica- 
tion IBM System/360 Operating System: Supervisor 
and Data Management Macro Instructions, Form C28- 
6647), such as CALL; or (2) coding assembler language 
branch instructions. 
In all cases, a save area must be provided that: 
© is aligned ona fullword boundary 
¢ is at least as large as the size specified in Tables 
14, 15, and 16, but preferably the standard 18 
words to ensure future compatibility 
° has its address in general register 13 at the time 
of the call macro instruction or branch 
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“Nore: Forrnan subprograms use certain floating-point 


registers (see Table 14), but do not save and restore 
original register contents. If the programmer wishes 
floating-point information retained, he must save it 
himself before calling the subprogram. 

If the called subprogram is one that uses FORTRAN 
input/output, error, or interruption routines (see Table 
17), the calling program must include the following 
two instructions before the branch is made: 

L ‘15, = V (IBCOM#) 

BAL, 14,64 (15) 
These instructions cause a branch into the 13com sub- 
program, which initializes return coding and prepares 
routines to handle interruptions. If this initialization is 
omitted, an interruption or error may cause abnormal 
termination. (After initialization, com returns to the 
intruction following the BAL. ) 


Nore: When these instructions are included, the occur- 
rence of a decimal divide exception within the assem- 
bler language program will cause the character B to 
appear in the program interruption message. 

When a branch instruction rather than a call macro 
instruction is used to invoke a subprogram, several 
additional conventions must be observed: 


e An arguinent (parameter) list must be assem- | 
bled on a fullword boundary. It consists of one 
4-byte address constant for each argument, with 
the last address constant containing a 1 in its 
high order bit. 

¢ The address of this argument list must be in 
general register 1. 

¢ The address of the entry name of the called 
subprogram must be in general register 15. . 

¢ The address of the point of return to the calling 
program must be in general register 14. 

The total requirements for an assembler language 
calling sequence are illustrated in Figure 1. 


Supplying Correct Paramefers 


Arguments must be of the proper type, length, and 
quantity, and, in certain cases, within a specified range 
for the subprogram called. 

For mathematical subprograms, this information can 
be found in Tables 2 through 6, mvTeceR*4 denotes a 
signed binary number four bytes long. REAL*4 and 
REAL*# are normalized floating-point numbers, 4 and 8 
bytes long, respectively. compLex*s and COMPLEX*16/ 
are complex numbers, 8 and 16 bytes long, respectively, 
whose left half contains the real part, and right half 





L 15, = V(IBCOM# ) 


BAL 14, 64(15) 
x e s 2 


LA 13, area 
LA L,arglist 
L 15,entry 
BALR 14,15 


NOP X’id” 


* + * * 


entry DC V (entry name ) 
or 
entry DC A (entry name) 
* * * * 
area DS xxF 


* e s s 


For one argument 
CNOP 
arglist DC X‘80’ 


DC AL3 (arg) 


For more than one argument 


CNOP 
arglist DC A (arg:) 
DC A (arg:) 
: 
° 
DC X‘80’ 


rere re > OTS ET SS OC SD LOCOCO 


These two statements are necessary if the called subprogram uses FORTRAN 


input/output, error, or interrupt routines (see Appendix B). 


General register 13 contains the address of the save area. 
General register 1 contains the address of the argument list. 
General register 15 contains the address of the subprogram. 


General register 14 contains the address of the point of return to the calling 
program, 


This statement is optional. The id represents an identification number. 
This number is supplied by the programmer and may be any hexadecimal 
integer less than 2" — 1. 


Note: In this case, the entry name must be defined by an EXTRN instruc- 
tion to obtain proper linkage. 


This statement defines the save area needed by the subprogram. The xx 
represents the minimum size of the save area required; however, the pro- 
grammer is advised to use a save area of 18 fullwords for all subprograms. 
(The minimum save area requirements are given in Tables 14 and 15 
for the mathematical subprograms and in Table 16 for the service sub- 
programs. ) 


Aligns the argument list at a fullword boundary. 
Places a 1 in the high order bit of the only argument. 


Contains the address of the argument. 


Aligns the argument list at a fullword boundary. 
Contains the address of the first argument. 


Contains the address of the second argument. 


Places a 1 in the high order bit of the last argument. 


Contains the address of the last argument. 


Figure 1. General Assembler Language Calling Sequence 


contains the imaginary part. Each part is a normalized 
floating-point number. Four-byte argument types must 
be aligned on fullword boundaries, and 8-byte and 16- 
byte types must be aligned on double-word boundaries. 

Argument information for nonmathematical sub- 
programs can be found in “Service Subprograms.” 

Error messages resulting from incorrect arguments 
are explained in “Appendix C: Interruption and Error 
Procedures.” 


Results 


Each mathematical subprogram returns a single answer 
' ‘of a type listed in Tables 2 through 6 (see “Function 
Value Type”). The mNTEGER answers are returned in 
general register 0, REAL answers are returned in float- 
ing-point register 0, and coMPLEX answers are returned 
in floating-point register 0 and 2. Result registers are 
listed by subprogram entry name in Table 14. 


The location and form of the service subroutine 
results can be determined from the discussion in 
“Service Subprograms.” 


Example 


To find the square root of the value in AMNT, the library 
square root subprogram (entry name sort) could be 
invoked, using the following statements: 


L 15, = V (IBCOM#) 
BAL. 14,64(15) 
LA " 13SAVE 
CALL SQRT,(AMNT),VL 
SAVE DS  18F 
AMNT DC "Baa? 
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Table 14. Assembler Information for the Explicitly Called 


Registers Used! 


(The vi_ operand in Cau indicates the end of the 
parameter list. ) 
Employing only assembler language instructions, the 


Mathematical Subprograms 


Save Area 


DMAXI, DMIN1 
DSQRT 

DTANH 

EXP 

ERF, ERFC 
IDINT, INT, IFIX 
MOD 

SORT 


Nw 
= ee 
a 


w 
i 
a 


returned in floating-point register 0 as a 4-byte floating- 
point number. 


Space Considerations 


Many of the mathematical subprograms require other, 
mathematical subprograms for their calculations. In 


Entry Name(s) (Fullwords ) | Result | Intermediate’ sequence would be: 
AINT ; 9 L 15,=V(IBCOM+#) 
ALGAMA, GAMMA 9 BAL 14.64(15) 
ALOG, ALOGIO0 7 ‘ 
AMANO, AMINO 6 . 
MANO, MINO 9 . 
AMANXI], AMINI 6. LA 13,.SAVE 
MAX}, MINI 9 LA 1,ARG 
AMOD. DMOD 9 2, 4.6 L IS,ENTRY 
ARCOS. ARSIN 10 2.4 BALR 14,15 
ATAN 5 0 2.4.6 . 
ATAN, ATAN2 7 0 2,-4.6 ‘ 
CABS 7 0,2 AG ' 
CCOS, CSIN 9 0,2 5 ENTRY ne V(SQRT) 
CDABS 7 0, 2 4.6 ‘ 
CDCOS, CDSIN ‘ 0. 2 4 . 
CDEXP 0,2 4.6 . 
CDLOG 0.2 AG SAVE DS 18F 
CDSQRT 0,2 4.6 . 
CEXP 0,2 46 . 
9 3 : 

COS. SIN ioe (ee ARG 1X: Lee 
COSIT, SIN 0 2.4 DC a 
COTAN, TAN 0 2.4 
CSQRT 0. 4.6 : 
DARCOS. DARSIN 0 2,4 . oer 
DATAN 5 0 2 4.6 AMNT DC Da 
DATAN, DATAN2 7 0 2.4.6 : 
DCOS. DSIN 7 0 2 4 Note that, in both cases, a branch to cox is pro- 
DCOSH, DSINH | 8 0 2.4 vided, a save area is set up, and AMNT meets argument 
poo aks DIAS : 2 Stc8 specifications by being a four-byte normalized floating- 
DERF, DERFC- il 0 2.4.6 F a ; 
DEXP 7 0 9 point number, = 0, aligned on a fullword boundary. 
DGAMMA, DLGAMA ll 0 2.4.6 In both cases, the answer (12, in this instance) is 
PDLOG, DLOGI10 ¢ 0 2.4.6 

0 

0 

0 

0 

0 

0 

0 

0 

0 


TANH 


‘Floating-point; asterisk indicates general. 





Table 15. Assembler Information for the Implicitly Called 


Mathematical Subprograms 


Save Area 


Registers Used! 


addition, most of the subprograms use the input/ 
output, error, and interruption library subroutines. 
(This interdependence is outlined in “Appendix B: 
Storage Estimates.”) Thus, although the programmer 
may request just one FORTRAN subprogram, the require- 
ments of that subprogram may make the resultant load 
module quite large. The sert routine, for example, takes 
only 344 bytes of storage itself, but requires other sub- 



















Entry Name(s) ( Fullwords ) routines that increase the load module size by approxi- 
CDMPY#. CDDVD# 0,2 4.6 mately 20,000 bytes. 
CDVD=, CMPY# 4,6 
TV p 
sae Table 16. Assembler Information for the Service Subprograms 
SAREE Entry Name(s) Save Area (Fullwords ) 
FDXPD2 DUMP, PDUMP 18 
FCDNIZ DVCHK 10 
FCXPIS EXIT 5 
OVERFL 10 


‘Floating-point; asterisk indicates general. 


















SLITE, SLITET 10 


Appendix B contains decimal storage estimates (in 
bytes) for the library subprograms. The estimate given 
does not include any additional mathematical subpro- 
grams or FORTRAN routines that the subprograms use 
during execution. The entry names of any additional 


Table 17, Mathematical Subprogram Storage Estimates 


Entry Name(s) 


AINT 

ALGAMA, GAMMA 
ALOG, ALOGIO 
AMAXO, AMINO, MAXO, MINO 
AMANI, AMIN], MAXI, MINI 
AMOD, DMOD 
ARCOS, ARSIN 
ATAN 

ATAN, ATAN2 
CABS 

CCOS, CSIN 
CDABS 

CDCOS, CDSIN 
CDDVD#, CDMPY# 
CDEXP 

CDLOG 

CDSORT 

CDVD#, CMPY# 
CEXP 

CLOG 

COS, SIN 

COSH, SINH 
COTAN, TAN 
CSQRT 

DARCOS, DARSIN 
DATAN 

DATAN, DATAN2 
DCOS, DSIN 
DCOSH, DSINH 
DCOTAN, DTAN 
DERF, DERFC 
DEXP 

DGAMMA, DLGAMA 
DLOG, DLOG10 
DMAX1, DMINI 
DSORT 

DTANH 

EXP 

ERF, ERFC 
FCDXPI# 

FCXPI# 

FDXPD# 

FDXPI# 

FRXPR# 

FIXPI# 

FRXPI# 

IDINT, INT 

MOD 

SQRT 

TANH 





Decimal Estimates 


80 
845 
464 
224 
224 
120 
496 
200 
488 
192 
760 
200 
$32 
240 

. G24 
488 
328 
216 
592 
464 
504 
50-4 
648 
312 
648 
312 
648 
696 
592 


~ 


é 
808 
704 

1056 
538 

120 
352 
304 
424 
520 
560 
536 

- 464 
368 
432 
365 
360 
136 

56 
344 
236 


Appendix B: Storage Estimates 


mathetical library subprograms used are given in 


Table 17. Tables 17 and IS also indicate which mathe- 


matical and service subprograms require FORTRAN 
routines for input/output. interruption, and error pro- 


cedures. 


Additional Mathematical 
Subprograms Used 


ALOG, ENP 


SORT 


SORT 

EXP. SINCOS 
DSORT 

DEXP, DSIN DCOS 
DEXP, DSIN ‘IDCOS 
DLOG, DATAN2 
DSORT 


EXP, SIX/COS 
ALOG, ATAN2 


EXP 


SORT 
DSORT 


DEXP 
DEX? 
DLOG, DENP 


DEXP 

ENP 

CDMPY= ‘CDDVD* 
CDVD*/CMPY= 
DEXP, DLOG 


EXP, ALOG 


EXP 








1/0, Error 
& Interrupt 
Routines 
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The programmer must add the estimates for all [This routine refers to a table (micuaTBL) for informa- 
subprograms and routines needed to determine the tion about the input/output devices used during 
amount of storage required. If the programmer has not —_—execution.] The m1cFCoME routine performs interrup- 
made allowances for the storage required by any of tion and error procedures for Basic ForTRAN Iv library 
these additional routines (see Table 19), the amount = subprograms; the mHCFCOMH, IHCFCVTH, IHCFINTH, 
of available storage may be exceeded and execution —_‘tHCTRCH, and 1HCUOPr routines perform the procedures 
cannot begin. for FORTRAN Iv library subprograms. If a system con- 
| tains both basic and full rorrraN rv compilers, the 
IHCFCOMH-IHCFCVTH routines are used. Tables 17 and 
System/360 Operating System 18 indicate which library subprograms require these 
The names of execution-time routines sometimes vary _ Cxecution-time routines. 
according to whether or not the Extended Error 


Handling facility is in effect. In the following discus- Table 18. Service Subprogram Storage Estimates 

sion, the names that are used are those when the facility 1/0, Error, 

has not been specified. Table 19 presents a cross listing & Interrupt 

of names for both circumstances. A full discussion of Entry Name(s) | Decimal Estimates Routines 

extended error handling is given in the FORTRAN IV DUMP/PDUMP Yes 

(G and H) Programmer's Guide listed in the Preface. Ds bas 
The IHCFIOSH routine performs input/output pro- OVERFL Yes 


cedures for both Basic FORTRAN IV and FORTRAN IV. SLITE/SLITET Yes 





Table 19. OS Execution-Time Routines Storage Estimates 


Decimal Estimate Used by 





Extended Error Handling Facility 
































IHCADJST 1,156 FORTRAN IV Same 
IHCDBUG 2,152 FORTRAN IV Same 
IHCFCOME Basic FORTRAN IV Not applicable 
IHCFCOMH FORTRAN IV IHCECOMH 
IHCCOMH2 FORTRAN IV Same 





IHCDIOSE Both IHCEDIOS 






(See Note 1) 
3,744 
+ IHCUATBL 













IHCFIOSH Both IHCEFIOS 


























































(See Notes 2 and 3) (See Notes 2 and 3) 
IHCFINTH FORTRAN IV IHCEFNTH 1,368 
LHCFCVTH FORTRAN IV Same Same 
IHCCGOTO Basic FORTRAN IV Not applicable 
IHCIBERH FORTRAN IV Same 
IHCIBERR Basic FORTRAN IV Not applicable 
» IHCNAMEL FORTRAN IV Same 
IE.CSTAE FORTRAN IV Same 
JHCTRCH FORTRAN IV IHCETRCH 
IHCUOPT FORTRAN IV Same 












(See Note 4) 
1,552 
824 




















Not applicable 
Not applicable 


FORTRAN IV 
FORTRAN IV 


IHCERRM 
IHCFOPT 






NOTEs: 


. This module also acquires dynamic storage. Its amount, in bytes, may be computed by the formula 184 + buffer size. 
Each buffer value should be 800, except for a card read punch which is 80 and a printer which is 133. FORTRAN utilizes 
double buffering. 
2. This module also acquires dynamic storage. Its amount, in bytes, may be computed by the formula 128 + buffer size. 
Buffer lengths are listed in Note 1, 
. The number of bytes in table IHCUATBL may be computed by the formula 12n + 8, where n is the number of data set 
reference numbers requested during system generation. 


. The number of additional entries supplicd in the Option Table during system generation is represented by n. 


In addition, several other execution-time routines 
may be needed to resolve external references in 4 
FORTRAN Iv object module. 7 
1. If a source module specifies direct-access input/ 

output operations, the compiler generates a call to 

the IHCDIOSE routine. . 

2. At the point that errors are encountered during 
compilation, the compiler generates a call to an 
error routine (mHiciBERR for Basic FORTRAN IV and 
IHCIBERH for FORTRAN Iv). If execution of the load 
module is attempted, the error routine is called, 
a thessage is issued, and the execution is terminated. — 


. If a Basic FORTRAN 1V source module contains a com- 


puted co To, the compiler generates a call to the 
IHCCGOTO routine. 


. If a FORTRAN Iv source module contains any input/ 


output operations that refer to a NAMELIST name, 
compiler generates a call to the 1HCNAMEL routine. 


. If a FORTRAN IV source module uses. the debug 


facility, the compiler generates a call to the nicbBUG 
routine. 


. If boundary alignment was specified during system 


generation, the 1HCcapysT routine will be loaded if a 
boundary-alignment error occurs. 
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Appendix C: Interruption and Error Procedures 


This appendix contains brief explanations of the pro- 
gram interruption and error procedures used by the 
FORTRAN library. A full description of program inter- 
rupts is given in the publication IBM System/360 
Principles of Operation, Form A22-6821. The program- 
mer's guides listed in the Preface give detailed informa- 
tion about error processing and message formats. 


Interruption Procedures 


The Fortran library processes those interrupts that are 
described below; all others are handled directly by 
the system Supervisor: 

1, When an interrupt occurs, indicators are set to 
record exponent overflow, underflow, fixed-point, 
floating-point or decimal divide exceptions. These 
indicators can be interrogated dynamically by the 
subprograms described in the chapter, “Service 
Subprograms.” 


2. A message is printed on the object error unit 


when each interrupt occurs. The old program 
status word (psw) printed in the message indi- 
cates the cause of each interrupt. 

3. Result registers are changed when exponent over- 
flow or exponent underflow (codes C and D) 
occur. (For a description of the format of floating- 
point numbers, see the publication IBM System/ 
360 Principles of Operation, Form A22-6821.) 
Result registers are also set when a floating-point 
instruction is referenced by an assembler language 
execute (ExEC) instruction. 

4. Condition codes set by floating-point addition or 
substraction instructions are altered for exponent 
underflow (code D). 


5. After the foregoing services are performed, execu- 
tion of the program continues from the instruction 
following the one that caused the interrupt. 


Error Procedures 


During execution, the mathematical subprograms as- 
sume that the argument(s) is the correct type. No 
checking is done for erroneous arguments (i.e., the 
wrong type, invalid characters, the wrong length, etc. ); 
therefore, 2 computation performed with an erroneous 
argument has an unpredictable result. However, the 
nature of some mathematical functions requires that 
the input be within a certain range. For example, the 
square root of a negative number is not permitted. 
If the argument is not within the valid range given in 
Tables 2 through 6, an error message is written on the 
object error unit data set defined by the installation 
during system generation. The execution of this load 
module or phase is terminated and control is returned 
to the operating system. However, execution can con- 
tinue, (in FORTRAN Iv (os) only), if extended error 
handling was selected during system generation. This 
facility provides for standard corrective action by the 
user. For a full description of extended error message 
handling, see the publication IBM System/360 FOR- 
TRAN IV (G and H) Programmer's Guide, Form 
C28-6817. 

Table 20 lists the error messages in numeric order, 
gives the entry name(s) associated with that number 
and an explanation of the error. The programmer's 
guides listed in the Preface of this publication show 
the text and format of the message as it appears in the 
error listing. In the following explanations, x repre- 
sents the argument supplied by the programmer. 


Table 20. Error Messages 


[ere [eerie 


SLITE, 
SLITET 


FIXPI# 



























An invalid sense-light number was detected in the argument list of a call to the SLITE or 
SLITET subprogram. 












In the FIXPI# subprogram, a base number of zero and an exponent less than or equal to zero 
is an error. 





FRXPI# In the FRXPI# a ee a base number of zero and an exponent less than or equal to zero 


is in error. 





FDXPI# In the FDXPI# subprogram, a base number of zero and an exponent less than or equal to zero 


is an error. 





FRXPR# In the FRXPR# subprogram, a base number of zero and an exponent less than or equal to zero 


is an error. 


FDXPD# 





In the FDXPD# subprogram, a base number of zero and an exponent less than or equal to zero 
is an error. 


FCXPI# 





In the FCXPI# subprogram, a base number of zero and an exponent less than or equal to zero 
is an error. 





FCDXI# In the FCDXI# subprogram, a base number of zero and an exponent less than or equal to zero 


is an error. 
SQRT 
EXP 


ALOG, 
ALOGI0 


In the SORT subprogram, the value of the argument is less than zero. 





In the EXP subprogram, the value of the argument is greater than 174.673. 












In the ALOG/ALOG10 subprogram, the value of the argument is less than or equal to-zero. 
Because this subprogram is called by an exponential subprogram, this message also indicates 
that an attempt has been made to raise a negative base to a real power. 









SIN, 
Cos 


ATAN2 


SINH, 
COSH 


ARSIN, 
ARCOS 


TAN, 
COTAN 


TAN, 
COTAN 


DSQRT 
DEXP 


DLOG, 
DLOG1O 


In the SIN/COS subprogram, the absolute value of an argument is greater than or ‘equal 
to gigs 





When the entry name ATAN2 is used, the value of both arguments cannot be zero. 









In the sppeeae SINH/COSH, the value of the argument is greater than or eauet to 
175.366 













In the ARSIN/ARCOS subprogram, the absolute value of the argument is greater than one. 






In the TAN/COTAN subprogram, the absolute value of the argument is greater than or equal 
to 238 x, 













In the TAN/COTAN subprogram, the value of the argument is too close to one of the 
singularities (+ /2, + 3~/2,... for the tangent; + wr, + 2x, ... for the cotangent). 














In the DSQRT subprogram, the value of the argument is less than zero. 





In the DEXP subprogram, the value of the argument is greater than 174.673. 










In the DLOG/DLOG10 subprogram, the value of the argument is less than or equal to zero. 
Because the subprogram is called by an exponential subprogram, this message also indicates 
that an attempt has been made to raise a negative base to a real power. 






DSIN, 
DCOS 


DATAN2 


DSINH, 
DCOSH 


DARSIN, 
DARCOS 





In the DSIN/DCOS subprogram, the absolute value of the argument is greater than or equal 
to 25° x, 





When the entry name DATAN2 is used, the value of both arguments cannot be zero. 









In the DSINH/DCOSH subprogram, the absolute value of the argument is greater than or 
equal to 175.366. 












In the DARSIN/DARCOS subprogram, the absolute value of the argument is greater than one. 
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Table 20. Error Messages (Continued ) 


Bomy Nene( 


DTAN, In the DTAN/DCOTAN subprogram, the absolute value of the argument is oe, than or 
DCOTAN equal to 2°° x 


DTAN, In the DTAN/DCOTAN subprogram, the value of the argument is too close to one of the 
DCOTAN singularities (+ w/2, + 3x/2,... for the tangent; + +, + 2x,... for the cotangent). 


CEXP In the CEXP subprogram, the absolute value of the real part of the argument is greater 
than 174.673. 


CEXP In the CEXP PLOT, the absolute valne of the imaginary bart of the argument is greater 
than or equal to 278 x 


CLOG In the CLOG subprogram, the value of both the real and imaginary parts of the argument 
cannot be zero. 


CSIN, In the CSIN/CCOS a ae the absolute value of the real part of the argument is greater 
CCOs than or equal to 218 x 


CSIN, In the CSIN/CCOS subprogram, the absolute value of the imaginary part of the argument is 
Ccos greater than 174.673. 


CDEXP In the CDEXP subprogram, the value of the real part of the argument is greater than 174.673. 


CDEXP In the CDEXP subprogram, the absolute value of the imaginary part of the argument is 
greater than or equal to 2°° x, 


CDLOG In the CDLOG subprogram, the value of both the real and imaginary parts of the argument 
cannot be zero. 


CDSIN, In the CDSIN/CDCOS subprogram, the absolute value of the real part of the argument is 
CDCOS greater than or equal to 2°° x. 


CDSIN, In the CDSIN/CDCOS subprogram, the absolute value of the imaginary part of the argument 
CDCOS is greater than 174.673. 


GAMMA In the GAMMA subprogram, the value of the argument is outside. the valid range of 
2 < X < 57.5744. 


ALGAMA In the ALGAMA subprogram, the value of the argument is outside the valid range of 
0< X < 4.2937 > 10”. 


DGAMMA In the DGAMMA subprogram, the value of the argument is outside the valid range of 
’ 2 < X < 57.5744. 


DLGAMA In the DLGAMA subprogram, the value of the argument is outside the valid range of 
0<X < 4.2937 10", 
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Table 21. Mathematical Subprogram Module Names 


@ Table 21. Mathematical Subprogram Module Names 
(Part 1 of 2) 


( Part 2 of 2) 


Module Name 





Module Name 









IHCCSABS 
JHCCSSCN 
IHCCLABS 
THCCLSCN 
IHCCLAS 

IHCCLEXP 















DTANH 


EXP 
ERF 
ERFC 
















FCDXI# 





IHCLTANH 


IHCSEXP 
IHCSERF 
IHCSERF 


IHCFCDAI 

































FORTRAN IV FORTRAN IV . 

Entry Name Entry Name Basic FORTRAN IV 
—_———-- IHCFAINT DERF IHCLERF | —----- 
THCSGAMA DERFC IWCLERF | -——---— 
IHCSLOG ITHCSLOG DEXP IHCLEXP 
IHCSLOG IHCSLOG DGAMMA IHCLGAMA | = -----— 
IHCFMAXI IHCFMAXI DLGAMA IHCLCAMA 100 3=<<=<s 
IHCFMAXR IHCFMAXR © DLOG IHCLLOG IHCLLOG 
IHCFMAXI IHCFMAXI DLOG10 IHCLLOG IHCLLOG 
IHCFMAXR IHCFMAXR DMAXI IHCFMAXD IHCFMAXD 

AMOD MIN} IHCFMAXD IHCFMAXD 
IHCSASCN DMOD So IHCFMODR 
IHCSASCN ISIN IHCLSCN IHCLSCN 
THCSATN2 DSINH IHCLSCNH | -—-—---— 
IHCSATN2 DSQRT IHICLSQRT 
DTAN IHCLTNCT | —-—---- 


IHCLTANH 
IHCSEXP 









IHCCLLOG FCXPI# TIHCFCXPL | £——-———— 
IHCCLAS FDXPD# IHCFDXPD IHCFDXPD 
IHCCLSCN FDXPI# IHCFDXPI IHCFDXPI 
IHCCLSQT FIXPI# IHCFIXPI IHCFIXPI 
IHCCSAS FRXPI# IHCFRXPI IHCFRXPI 
IHCCSEXP FRXPR# IHCFRXPR IHCFRXPR 


IHCCSLOG 
IHCCSAS 
IHCSSCN 
IHCSSCNH 
IHCSTNCT 


IHCSSCN 


ee ee oe eer eee Som 










IDINT 
INT 













GAMMA 


IHCSGAMA 








IHCIFIX 
THCIFIX 


oo MAXO IHCFMAXI IHCFMAXI 

IHCLSCN IHCLSCN MAXI IHCFMAXR IHCFMAXR 
IHCCSSQT = MINO IHCFMAXI IHCFMAXI 
MINI IHCFMAXR IHCFMAXR 

IHCLASCN MOD -——— IHCFMODI 


IHCLASCN 
THCLATN2 
IHCLATN2 
THCLSCN 

IHCLSCNH 
IHCLTNCT 


IHCLATAN 


IHCLSCN 












SIN 
SINH 
SQRT 


TAN 
TANH 






IHCSSQRT 












IHCSSCN 
IHCSSCNH 


IHCSTNCT 
IHCSTANH 
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IHCSSCN 












IHCSTANH 














Appendix E: Sample Storage Printouts 


A sample printout is given below for each dump 
format that can be specified for the storage dump 
subprogram. The printouts are given in the following 
order: hexadecimal, LOGICAL*1, LOGICAL*4, INTEGER*®, 
INTEGER*4, REAL*4, REAL*8, COMPLEX*8, COMPLEX*16, 
and literal (see Figure 2). Note that the headings on the 
printouts are not generated by the system, but were 


CALL POUMP WITH HEXADECIMAL FORMAY SPECIFIED 


obtained by using rorMaT statements. The number 
printed at the left of each output line is the storage loca- 
tion (in hexadecimal ) of the first data item tabulated. 

The output of the storage dump subprogram (for 
both entry names, pump and poump) is placed on the 
object error unit data set defined by the installation 
during system generation. 


OOASED SOSFSELD 00000000 SESFSEIO 10000000 %2100000 

oo6ncs %2800000 00000000 00000000 00000000 00000000 00000000 96000000 ococne0e enndened oefeccnD coneonoo Bsoecoeoje 
oo60F8 €0000000 00090000 %1200000 S1566666 do00000c &1200000 oat 

CALL POUMP WITH LOCICAL®L FORMAT SPECIFIED 

OOGELE % fF 


Ne SY SD SR crt aS 


CALL PCUMP WITH LOGICAL®R FORMAT. SPECIFIED 
006£10 fe Tf 


Ae eS a ST RES 


CALL POUMP WITH INTEGER®2? PORMAT SPECIFIED 


COGEI8 10 
OOGERA 100 
oCGerc 10 


CALL POUMP WITH ENTECEROS FORMAT SPECIFIED 
006220 


1 2 5 & s 

OObESs 13 : 12 ee 
CALL POUMP WITH REAL®® FORMAT SPECIFIED 
6CGEOO O.20000000E 01 0.939999968 01 =e 
CALL PDOUMP WITH REAL®S PORMAT SPECIFIED — : 
ooboCs 0.17539999999999999D 03 = 
CALL POUMP WITH COMPLEX®S FORMAT SPECIFIEO 

_ 06000 €3.¢000000,%, 00000005 


CALL PCUMP WITH COMPLEX®16 FORMAT SPECIFIED 


oceoee €0.9999999999999990, 0, 9999998999999990) 
CALL PDUMP WITH LITERAL FORMAT SPECIFIED 
OoGEsC THIS ARRAY CONTAINS ALPHAMERIC DATA 


Figure 2. Sample Storage Printouts 
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(%,0000000, 8.0000000) 


ee Nee ce mcnaeriee 


(~0.9999999999999990, -0,9999999999999990) 


absolute error ........0.......00.. wiht anicnates, te toe 18, 43-47 
absolute value subprograms ..................-.. .6, 10, 20 
accuracy figures for mathematical subprograms... .... 43-47 
AINT subprogram 
oe sd gite CUNT Se S TE SR a aoe ew ehh Rd 55 
sana iey ace inate te van Ont Darth tela cee eg indy Sars ll 
ALGAMA/ GAMMaA subprogram 
ACCUTACY 3a e 5 cc Pees eb REE dead we 44, 46 
SIGOLNIM foi 95 6 8455 CRRA AS ee ee SNS 27-28 
effect of an argument error ....... ............5. 28 
CPOE MOSSARES iin ce cid Hae yaa Sian eee he wa 60 
SIZE. Sena hor ain de hie Be enee oi Gaon eho 55 
MEMNSS “esd 5 ars wduniulrpeds Sarees pate ese eee ees 48, 50 
NGO vein Cale Sees et wet yi Rea nies Bee edies 10 
algorithms for mathematical subprograms ............ 17-42 
ALOG/ALOG10 subprogram 
BCCULOCY 2 fare PSN Bi ode OL Alt ea OR aie 44 
BUDOMIRON «tance Se AO ea oct etmek ee ee See . gil 
effect of an argument error................... 31 
@ELOYK MESSARE: iors eke By Rea URN aes 59 
SIZE olathe Beale ES Hee EOS Br ets 55 
ella sc Dive uh SOs RRS deres tse ta Dhacda tite ity aes Areca ae ana 48 
Srp deus Ace tstedean tee otk otis Wiegand Dae a oet Rae Deak te eh 7 
AMAXO/ AMINO/MAX0/MINO subprogram 
— 5 Sea Seal poate an ete dene ee ee Aad oh erntagt ane areie aoe 55 
ote eRe gh ase te eee walirde tah hee Sed artes 28 a Bk alte hab aves 11 
AMAXI/AMIN1/ MAXI/MINI subprogram 
— une ee ae ee deb oy siete Les wedges er 55 
Macro ia ttea ise alate mate Scares bak ee ee eae 10-11 
Assnieas National Standards Institute (ANSI) ...... 
AMOD/DMOD subprogram 
WIZE 2 ive hed aas Gas Sine Gad oy Seu yee ee 55 
UGG 68 is Faw Ao cede. side Gee bike de wn a SN as ll 
arccosine subprograms ...............0 .cceeeeee 6, 8, 20-21 
ARCOS/ARSIN subprogram 
QUCUIACY 5 35 Gbawe ee ie lass cr Rea ea PO as 44 
BIQORBAE ib. heb hk ORE Srrere le Dad pee are ROO 20-21 
effect of an argument error .................00.00: 2) 
CMVOF MESEARE i's Si ges Raise we ee ee ev age Vea 59 
WIZE chon atta Palen bee eae Ae ORT A ees 55 
Hine eric ace dics cued busins ta taeee Pisin nal e 48 
MSG ecu Bh dase hea Ne ee CEO Res See Dees 8 
arcsine subprograms ...............0 6.005. :6, 8, 20-21, 55 
arctangent subprograms ..................... 6, 8, 22-24, 55 
APRUMEHS one owt eh a sew a et en eE ewe 8, 53-54 
ARSIN/ARCOS (see ARCOS/ARSIN ) 
assembler language calling sequence ................ 52-54 
assembler requirements .................-000 0002s 52-54 
.ATAN (Basic FORTRAN IV) subprogram 
ROCUINCY cities hie SL eh ada eh 6 2a eplenee ee 44 
RIQOTEIN o6cihie ea ole cice Mewes sites Gp Pee 22 
effect of an argument error ...................... 22 
37 = ae oO Pa a ee eer 55 
a Srey, Reales SW te oat uh eh eesti aed tard el chasse 48 
sites tte uc thcetamstatta Ciubld atx cticay etre tsanara aban atts coats cand en av 8 
ATAN/ATAN2 (FORTRAN IV) subprogram 
ACCUIACY: Foci ote eae sabe ae tee ey east 44 
ISOPEN 55555 Foch Bie pede she ued ach ade daraigiutecinis 22-23 
effect of an argument error ...................... 23 
@IYOr MESSARG «iis coe das Rd daw 6 oN akg ox 59 
SING 5 bik Ors mate es Hel BAS EGON Whee A REE 55 
MUMS 2 5h atecee wth pans ou eked be ewe eae anereee 48 
MRO Fests ie cise at Sh a a er en DE 8 
CABS subprogram 
BCCUIRCY Ss ccs Bast ene san di ee aad atoen peas 44 
elgoritlim: occas sls ee sarah esd cree nabee ens 20 


Index 


effect of an argument error... 06... ewe ee 20 
SIZE ide. thn TA ee, eh eg ity ected Sane: 55 
CHINES: Af oe. B.S See gescsnhoenee ta eres 48 
MISE aches 3 seq Rc aiersen Se ue Mee nea: * ER: Ga 10 
calling sequence in assembler Janguage. .......... . §2-53 
calling FORTRAN subprograms 
EXDIGHY 3 ou vag Gea GLERY He TENE bee eiod 7-11 
implicitly ...... REL Meets Gah 2 Shien acd asked. ph Ge a BE -... 12-13 
in assembler language ........-.......0.00 000005 52-53 
CALL macro instruction 2. ........220 00.5. ce aeee 52 
CCOS/CSIN subproerem 
accuracy 2.0.2... bah tadec unis pleat chien canes 44, 45 
AIPOHNM iG. Ss. ° 23. of Mrtagins QeSeey o Fas 35 
effect of an argument error .......00.0 2. -0- ee eee. 35 
OYTOR WESSALES. oe es bacea Pasi ae aa 60 
BIO teens. 2 hare aenvadoveoneaacwante 55 
HUDINGS: aceneh, ge dele . BPRS ete yes eG 48 
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accuracy =... ee er re rn ee 44 
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effect of an argume nt ¢ Pett) eee CRE ere re 20 
SIZ00 3 Scie Basen cen GPG Ahan at Ee 55 
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CDCOS/CDSIN subprogram 
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size... ......- See er ne nets a eee 55 
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CDDVD#. / CDMPY # subprogram 
restilt OF SE: i be hae ee ae ew tees ae 12 
WPS cs vos, woe Medd ARE PMOL U ES Rela hy aeeS ees 55 
use 
CDEXP subprogram 
RCCUPACY: 8G a AeA: Sewioe es HEE ERY 44 
HI GOTICHIN oe eS Se, -Eearaie tee a eas ele 27 
effect of an argument error .......0 00... eee ee 27 
Error MESSAKES ©... eee eee 60 
CL) ee eee ee eer cre are eres 55 
timings: © sk. Sa weed Se wowed bese eeses 48 
Sts A bine via Aa en ed ee eles! ene ea ee Data of 
CDLOG subprogram 
ACCITACY 2c es Gabi Odes ede ble sea sR ey 64 5'3 44 
algorithm ..........-... 0... het sg apie a pa tates 32-33 
effect of an argument HOR. ot io oe nee se es 33 
error message .0 0 6.6 eee eens 60 
SIZE Meiers cee. oh bee MS Se BESO Gale ark ibaned eke 55 
HININGS Arg oweua iene ead tnlseeoong ad eta 48 
NSO oie iy AE aee phatase Shoe ey Vik aes 7 
CDMPY#/CDDVD# (see CDDVD#/CDMPY#) 
CDSIN/CDCOS (see CDCOS/CDSIN ) 
CDSQRT subprogram 
MOCUTACY ices Sats Coes Poteet Ges Cea Cae 44 
{ @IGOMENON 2250. 4s.n gee aKeaw ea ha ea Oeen wie ee earnest! 37-38 
effect of an argument error .......0 0.62 -- 00sec eee 38 
$708 i oe ata eancthe eae ae has Bae a deg ES 55 
GintngS hed enti e ee eis ested te sa toys et one -48 
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CEXP subprogram 
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timings .. ....... ePtcne ohtehreas draieetetan a Sean 48 
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CLOG subprogram 
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RIZG cme Ras esd St pene leg np anh eM eg 55 
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CMPY#/CDVD# (see panes A paae 
common logarithm subprograms.... . . 6, 7, 17, 27-29, 31-33, 55 
complemented error function subprogram ...... 6, 10, 24-26, 55 
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COS/SIN subprogram 
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effect of an argument error, ........ ...........-. 34 
CIIOL MOSSARE™ nn yes iin th oir ek pes eee Saw Meek 59 
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effect of an argument error .................... 29 
CITOF MNUSEARE: foe deck tala .waterees a yce tie 59 
SILOM cca sete, sfiakinw en “cbatesad Gna ten 55 
RUPEES 25. Se wi coeds Wig Se, Ma meets ic en eden ane om 48, 50 
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cosine subprograms .... ............. .... 6, 8-9, 33-36, 55 
COTAN/TAN subprogram 
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SONI EB aco cles serps te Gio c ae Wh Sse donee 38-39 
effect of an argument error ....... se wipe ssenta te 39 
error MesSageS ... 60... cece ees 59 
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TUN GS a eg eh acng ocd we EG eed bameles ctec on eee 48, 50 
MSE! 25. Gabe tutes DORM eRe Rees ela OEM a 
cotangent subprograms ...................... 6, 9, 38-40, 55 
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CSQRT subprogram 
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SIZE? ito Setucd Ng a doy EN ey Beane oe rete a anwlala 55 
timings. ............., rte eco tnde, varied Gaus 1% ean 48 
MISO ig seh ah.ar itech te dot Gov era dat sy Seok ig an eget hes, 7 
DARCOS/DARSIN subprogram 
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BYRON 5 ulna dein aa hat Mohs Oe alana ney heats 34-35 
effect of an argument error....................05. 35 
EMrOr MESSAGe .... 1... ee eee ee eee eens 59 
SILC oii h 2 ocr ach dim mtr et ane eae task Ue oles copra tits 55 
CINTIB ES 3-<* \igeik. Gh sired be Boers e's asks Ole wale 48, 49 
MISC Sscg Sxl Sota a eerie ee hha ia ore nee a Led 
DARSIN/DARCOS (see DARCOS/DARSIN) 
DATAN (Basic FORTRAN IV) subprogram 
COCO iy taacenc ds Uh ences ERS ei le ae Werne en ete Mok 45 
BIPOHEANY, cscs cds wis eck aang Obed win enhack how aceguees 23 
effect of an argument error ..................... 23 
SFO Lien des aa he Gace Ne eg LN eee aeRO Bhan 55 
MMINGS oof tts tars eta Sawa peele wake teed eas 49 
WSO ei Bey dee wai os Sn SE hk ed 8 


DATAN/DATAN2 (FORTRAN IV) subprogram 


ROOUWACY 22cm erin e bas ba te nea oO ee 45 
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effect of an argument error .. 2.0... eee 24 
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size.......... PERE REO Oe eRe Te eee - §5 
tonnes h S Sen cosine eyraiyGi arama bhsainise a eaten ee take on eae 49 
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DCOS/ DSIN subprogram 
MOUUINC «5.8 eas boune Daud heen pee Ee even NOs 46 
RIVOMENO Ges icy heed te See ees 34-35 
effect of an argument error 20.0... eee 35 
(FTOT ANESSALE 6 bo es on ae vduses eee eis eaeee es 59 
SIZE? like ne edee ky eter oe Danang a RE Ss een are ae 
oe , deat dic nepcgtin a Bratt elesdcstaey apgecks oer ane ar aki ee 49, 50 
DCOSH/ DSINH subprogram 
accuracy ....... ee es eee 46 
algorithm ............ Gai Monchiad team iy saat 29-30 
effect of an argument error ........0..00.0.0 0000 ee oe 30 
OFTOE WOSSARE: igo 55 RG EE REY LR eRe PERE 59 
S126 Foi Ph areata a aid OG ath whee Pee eae nee 55 
ae Aes Ske a OU ger an De ate rat lead Se 8 ahd ens 49, sa 
DCOTAN/ DTAN subprogram 
AOOUPROY 5 dedstcae daw ne one ah Olea car edna easton wee 45 
RIBOTINIG: elie claliv ae Gaxnants sawed ataguaaia sete 39-40 
effect of an argument error ............0. 0002 eee 40 
@TTOT MESSAGE 6. boeken ek pada ewe et od no ae dees 60 
SIZE bia ede nd Semen s dae seas Maaeeatins 55 
ae Bian A snp Bote Ote le Rasa eaa tart waa amie ates pee Nal 49, 50 
DERF/ DERFC subprogram 
NOCULACY cay gs 4cda a Sees ee tweens 45 
OIMOMM NN: 56s bah foe ee Sea Ee Ope eT eee 25-26 
effect of an argument error ............00000eeeee 28 
SG oe ee as et ce eea a sale det Sas 55 
timings .............. ‘stata ore Whaat eee ek ee ere 49 
MSG ceetic df PY EAS a ws clans ain aieg See Wn eo ary aah a saa 10 
DEXP subprogram ’ 
PECULACY o/c bse/p sy Saws SG NR SRS a Oe ee hase ee eae es 45 
BlGORENN 54 5 abe aire Saws Shoe baw eases 26-27 
effect of an argument error ........ Eee eon es 27 
ETTOT MESSABE . 0... ccc cece eee encase 59 
C1 se ae ai A eee OL LON ENCE Sk mea ear 55 
aad OE ee eT ry eee Tree ere ne 49 
DCAMMA/ DLGAMA subprogram 
HeCUCY 635 ite ls hana beams aa anes 45 
BIPOLIUNIN, d26:2 oi Soe arb aoe acne wad oaswataeaede 28-29 
effect of an argument error ................0.0055 28-29 
CIFOF MOSSAPE: ov iis wa ba Te Chad See ee oe wae es 60 
BIZ Saath SAM GAY ts ee thn re ata ein ate leeds 55 
ag Shee te i th Saul Wg ec aac ee eh yee ae 49 
a ot tinaiglie § aon ghia hetecgn @ aek gaa area aM wie acted 10 
divide-check service subprogram .............-.25: 14, 54, 56 
(see also DVCHK) . 
divide-check exception ...............ceeeeeeeeeees ll 
DLGAMA/DGAMMA (see DGAMMA/DLGAMA) 
DLOG/DLOG10 subprogram 
ACCUIGCY” <hr dn ney haenooew uae ated Nota ioaiaens 45-46 
ASOT bs 4 uve Ss i case eeuawandet nee: 32 
effect of an argument error .............-..000005 32 
@MVOr MeSSAZe ... 0... ee ccc cece ee te eens 59 
SO Se ce ek Awa aA EOE 55 
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DMAX1/ DMINI subprogram 
1 Ee ee TE ree weige ais 55 
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DMOD/ AMOD (see AMOD/DMOD) 
DSIN/DCOS (see DCOS/DSIN) 
DSINH/DCOSH (see DCOSH/DSINH) 


DSQRT subprogram 


ACCUTACY Ye ines Sate tyeny Se Ste eek? pa bte ot eee eae's 46 
RIGOHEAN eo ees Pe Eons ies 36-37 
effect of an argument error ..................0005 37 
CTTOF TNCSSAQE, 6 es soar a eae RG ee ses es 59 
BAD sean Sie and a Es tie has PEE RC A CO 55 
CMSs oso eco ak whe Se ea bien 2 sere ase e 4 50 
MSE - 3io dicho ot ds Peas MARRS he Se kee eee 7 
DTAN/DCOTAN (see DCOTAN/DTAN ) 
DTANH subprogram 
ACCUIACY oe so ED oss serarmeattc Sate 46 
algorithm .........0.0..0 00.00 ccc cee ec eee ee eee 30-31 
effect of an argument error ..............-....5.. 31 
BIZ in Hane tan eure nd. ncaanilomea Mean aee ew adads 55 
MIU RS 365 tae aides oa ead ome se re tana dd wires 50 
MSG. 2 icptirise pie hah hee Ls bo ielnte aaah a eee anwar 9 
DUMP/PDUMP service subprogram 
assembler language requirements............. .... 54 
format specifications ...............0.0..2.20.00ee 15 
OUROUE: sos Pic P talline hy day sla a ok wae ne Soneeree 62 
programming consideration .....................- 15-16 
sample printouts ...... tPasie nate eye uaseenasG 62 
BEL is She os os Pie wl amy ese Bal ld Stele Re eek sais 55 
MSG oe in hatha ct ata tse ee Marais Shits Eo iy geese ne 15-16 
DVCHK service subprogram 
assembler language requirements ................. 54 
6 Sead hd ddd OR Ds ne ee 56 
SC thee 6 tobe oe. ore ads tend wae ined 2 rete ula ace 14-15 
end of execution service subprogram ... .......... 15, 54, 56 
(see also EXIT) 
Ontry NaAMe scoo5 2 hk, oe Ved oe eh he wd eee 7 
ERF/ERFC subprogram 
DCCULACY sigs Cosette bia. caw gets Pw tnx aces 46 
SORE IN oro 9) ances wih Ate odes beste ak hee Sis eS ore ey 24-25 
effect of an argument error ...................... 25 
SIG ei Rit ok We eR Been nak te ba aes 55 
RMINGS 3.b5 os Bobet es 4 of olka Wen bn Bae yews 50 
MSG ie scrsch Gras Pale ae ena Wien MURS ata Berea iso 10 
error 
oie) (a a ea a a 18, 45 
INCSSARES | 025 oes og ahacbainess Se awe Reels aa Sed 59, 60 
optional service (see extended error handling facility) 
DIOCEQUIOS ¢ S5 biotch thei a al ee A et ae 58-60 
propagation .......... 0.0... ccc cee ec eee eens 43 
POIAUVG 55925 0h cB bul as ealewe eno Neaes ats 18, 43 
error function subprograms ................. 6, 10, 24-26, 55 
execution-time routines .............0.. 0.020000 cue 55-57 
EXIT service subprogram 
assembler language requirements ................. 54 
SIZ ro eadhn™ aie ccne bipth ois Palen dca Ma ae Ouest ee 56 
USE sa see seis ied fale: ©. ek Ge yA ed Sa ccotite we a RVaneneed bb alehen e aAON WES 15 
EXP subprogram 
GCOUIROY 5 cs Sows es ec eek ba See Nae doa ERS aan 46 
lori tang 5 ere 5 so 4 HA hte KE lee Seek ooo 26 
effect of an argument error ..................-... 26 
CITOF MESSABE 2.0... ccc ecco eee acco acectces 59 
ZO ice rea k ewes eyte pte elo shia Hee We eee ls 55 
MINE eet ete ne BW Sel ot He costa Bb a Mig do ed 50 
NG ink eaeinee slik cme mag nee mee ene dus 7 
explicitly called subprograms ...................... 5-11 
accuracy figures .........,..... 000. cc cece eceees 44-47 
MABE 3 soicnia wre a ay oa tehesiialce accra sted acide eo aoe 6 
performance StatisCS. fos nce caw ee’ Der twdeat ‘43-51 
BIE sche pty esas ai Uk Sh wee ited aes 55-57 
MEDICS hatch ad Sapa craut aaa eenekeuh aa hee 7-11 
timing estimates ............00..0..0 cc cece eee, 48-51 
use in assembler language ....................... 52-54 
use in FORTRAN ................. ccc ccc e cence 1-8 
exponential subprograms 
OXDUCIE edn 2 3c ehh PRA bade 6, 7, 26-27, 55, 59 
MMHG os ass hecae yun wlekeaeon se ee 12-13, 40-42, 55, 59 


exponentiation 
explicit (see EXP) 


implicit 
with complex base and integer exponent.... ..... 13 
with integer base and exponent . .............. 13 
with real base and exponent .................. 13 
with real base and integer exponent ...... ..... 13 
exponent overflow exception. ....... 000.22. 0.0005 14, 58 
(sce also OVERFL) 
intermediate: < acs “Saiyaikaw ddd cud ee deeds Rast 18-19 
POP 633 aad ie Ae: Use oe eee Olds 18-19, 58 
exponent underflow exception di Gaiea gua Dnebarnce leah oe 14, 58 
(see also OVERFI. ) 
intemi@diate: c.5 4 Deasksk Saleh aedeseerssat 18-19 
terminal 65.0, 348 ta carsee? ava snes Ieee 18-19, 58 
extended error handling facility .................. 14, 56, 58 
FCDNI# subprogram 
error message ...0 008 eee 59 
result of me@s8. 6 (AS aa Daeaawt ea eye ten ee eens 12 
Sie Ek ae | ek aha Oto Be bien aad 55 
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FDXPD# subprogram 
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MIZE ie ht ees <a oh bun cca bea ieeea oes 55 
MSG? oof a Me eas Unt eteseite Date aah oh toss Wn acedia eceeeee 12 
FDXPI# subprogram 
error message 22) << 22: dase Fey sci enae es eet C943 59 
result of use ....... = siviee ecu toe teh Mink ea Pact aes 12 
a S Sais eats taste ee hp tee ethane oR kinins eet ahs 55 
PORES Pere ee eer ee Pe ee eee Se 12 
F IXPI# subprogram 
CIOL -MESSARE 2 See 65 Se EAs ea 59 
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SIZO 5 sos ae nk a das ee eG ee Set oe ea ae 55 
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FRXPI# subprogram 
@ITOR -MeSSACE 65-5". 5 how Ra Ns eee eas 59 
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ta Seis oid wietics silos stb te Sb as gts esa rn tin nenas ae 55 
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FRXPR# subprogram 
CIOL MESURE 56.04. co Cees wa dee eS . 59 
pestilkOF USO). co seas! oe eis dee ieaew Seeeeas 12 
BIG ea ou ab Mie a tle aK EE abe smn EMA ase 55 
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FUNCTION statement ...0..0.... 0.0. cee eee 5 
function value ..... 0.0.0.0... 0c cece eee ee es 5, 7-12 
GAMMA/ALGAMA (see ALGAMA/GAMMA) 
gamma subprograms ...........-..5500.000% 6, 10, 27-29, 55 
hyperbolic cosine subprograms ....... ........ 6, 9, 29-30, 55 
hyperbolic sine subprograms . .............. 6, 9, 29-30, 55 
hyperbolic tangent subprograms .............. 6, 9, 30-31, 55 
TBCOMSE: 60553505 -. seaetaisai aperears SBS ne 52 
IDINT/IFIX/INT subprogram 
= ORE itaca ele PODS te Bed eatin Beat PERN gee hates 55 
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IFIX (see IDINT/IFIX/INT) 
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THOCLADS (iS ces fet Bo ae ea ois eR Ras 61 
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IHCCLAS (see CDODVD# ces ) 
IHCCLEXP 
(see also CDEXP) 
THCCEEOG: 2. 3iiiseser ike’ wikien % aria Heced 
(see also CDLOG) 
IHCCLSCN ......... ..... 
(see also CDCOS/CDSIN ) 
IHCCLSQRT . . 
(see also CDSQRT) 
IHCCOMH?2 ....... 
IHCCSABS ._.... 
(see also CABS) 
IHCCSAS (see CDVD#/CMPY# ) 
IHCCSEXP.... ... wee eee 
(see also CEXP) 
IHCCSLOG ..... 
(see also CLOG) 
IHCCSSCN 
(see also CCOS/CSIN ) 
IHCCSSQT .......... ..... 
(see also CSORT) 
IHCDBUG . 
IHCDIOSE 
IHCECOMH 
IHCEDIOS . 
IHCEFIOS 
THCEPNTH: . «fo Joce28o.eben kek) Sead taseaes 
IHCERRM ......... ... 
IHCETRCH__.. ......... 
IHCFAINT (see AINT) 
IHCFCDXI (see FCDXI# ) 
IHCFCOME.... 
IHCFCOMH ... 
IHCFCVTH ...... ........... 
IHCFCXPI (see FCXPI# ) 
IHCFDXPD (see FDXPD# ) 
IHCFDXPI (see FDXPI# ) 
THCFIFIX (see IDINT/IFIX/INT ) 
IHCFIOSH ............. 
THOPIN GE 3 elcx 2etuecndtwoancee 2b ate Seng cer 
IHCFIXPI (see FIXPI# ) 
IHCFMAXD (see DMAX1/DMIN1) 
THCEMAXI (see AMAX0/AMINO/MAXO/MINO) 
IHCFMAXR. (see AMAX1/AMINI/MAXI/MIN1 ) 
IHCFMODI (see MOD) 
IHCFMODR (see AMOD/DMOD) 
THOEOPT, oie cy 2a ok Hana Op werent 
IHCFRXPI (see FRXPI# ) 
IHCFRXPR (see nit dal 
IHCIBERH . ... 
IHCIBERR....... 
THOLASON tcl ee he go Boa es eaten 
(see also DARCOS/DARSIN ) 
MHCLATAN? sociiteb SG oa kates 2h Oa tek lo b4 
(see also DATAN-Basic: FORTRAN IV) 
PECL) © 52 Sosa idnevaces Ba press pete le Se knw lend 
(see ee DATAN/DATAN2-FORTRAN Nay) 
IHCLERF .................... 
(see feo DERF/DERFC) 
FHICLEAP eect cue cikiean ham etaee Poet whaeaela. deat 
(see ae DEXP) 
IHCLGAMA .................. .... 
(see alto DGAMMA/DLGAMA) 
WHCLLOG ............ 
(see alsa DLOG/DLOGIO) 
WICESON 5,655. is tak Satis 
(see also ao DSIN) . 
KHCLSCNH ............ 
(see ie eee DSINH) 
THCLSQRT .... 0. ccc ne ee 
(see also DSQRT) 
THOR TANG eos vcrena pales acd! Mion Sw ewtieae es 
(see also DTANH) 
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THOSASON:: jdeckcne car ad Sess VG RANE 3 Sew anes 61 
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THOSERP fcc or. bee ea eee sae eee 4 Bee 61 
(see also EXP) 
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PHCSEOG 5 5.5 biscaia 8 Peau. BESS WA OU tad 61 
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